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1. 


INTRODUCTION 


1.1  GGSS  BACKGROUND 

The  Air  Force  Geophysics  Laboratory  (AFGL)  and  Defense 
Mapping  Agency  (DMA)  are  pursuing  the  development  of  a  moving- 
base  Gravity  Gradiometer  Survey  System  (GGSS).  The  GGSS  will 
use  gravity  gradients  to  improve  greatly  our  knowledge  of  the 
earth's  gravity  field  and  to  increase  significantly  the  rate 
at  which  we  obtain  that  knowledge.  Plans  call  for  incorporation 
of  this  new  gravity  sensor  in  both  a  land  vehicle  and  a  fixed- 
wing  aircraft.  Surface  testing  is  expected  to  commence  in  the 
spring  of  1986;  airborne  testing  in  the  fall.  Following  com¬ 
pletion  of  flight  testing,  deployment  of  the  GGSS  for  transi¬ 
tional  test/survey  operations  will  occur  immediately. 

Updates  from  the  Global  Positioning  System  (GPS)  will 
be  used  in  the  aircraft  application  to  augment  the  positions 
which  are  derived  from  the  sensor  outputs  of  the  GGSS  inertial 
platform.  Since  there  is  a  potential  for  significant  misregis¬ 
tration  during  surface  testing,  updates  from  wayside  position 
markers  and/or  a  fifth  wheel  are  being  considered  by  Bell 
Aerospace.  An  overview  of  the  currently  configured  GGSS  is 
presented  in  Fig.  1.1-1.  This  overview  also  indicates  the 
two  GPS  back-up  modes  which  are  a  part  of  the  baseline  GGSS 
configuration.  The  first  back-up  mode  is  initiated  when  only 
three  satellites  are  visible.  In  this  case,  incorporation  of 
an  atomic  clock  allows  the  full  navigation  solution  (i.e., 
three  positions  and  time)  to  be  computed.  The  second  back-up 
mode,  which  involves  the  use  of  both  an  atomic  clock  and  alti¬ 
meter,  will  be  called  upon  when  only  two  GPS  satellites  are 
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Figure  1.1-1 


Baseline  GGSS/External 
Reference  Configuration 


1.2  STUDY  OBJECTIVES  AND  REPORT  OVERVIEW 

The  primary  objective  of  this  study  is  to  identify 
and  analyze  navigation  system  configurations  which  are  most 
appropriate  for  assuring  the  immediate  and  long-range  success 
of  the  surface  and  airborne  GGSS  applications .  To  accomplish 
this  objective,  alternative  gravity-  and  navigation-aided 
subsystems  are  investigated;  current  and  anticipated  accur¬ 
acies  appropriate  to  each  subsystem  are  determined;  and  the 
resulting  accuracy  and  other  benefits  which  would  accrue  from 
using  these  aides  with  the  GGSS  are  assesssed.  The  following 
intermediate  study  objectives  are  also  addressed: 
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•  Review  sensor  technology 

•  Formulate  appropriate  error  models 

•  Examine  relevant  mechanization  issues 

•  Calculate  resulting  benefits  (in  particu¬ 
lar,  accuracy  improvement) 

•  Identify  technical  factors  to  be  overcome 

•  Recommend  areas  which  merit  further 
research . 

Four  classes  of  survey  aides  are  evaluated  in  this 
study  in  terms  of  their  ability  to  provide  significant  accuracy 
improvement  compared  with  the  unaided  GGSS.  The  first  is  con¬ 
ventional  GPS  and  its  associated  back-up  modes  (Chapter  2). 
The  second  involves  radio  interferometry  u* ing  GPS  signals 
and  a  separate  GPS  receiver,  atomic  clock,  weather  sensor, 
and  data  recorder  at  a  convenient  fixed  location  (Chapter  3). 
The  third  is  a  gravimeter  which  provides  an  additional  (and 
unique)  gravity  measurement  (Chapter  4).  The  fourth  is  a 
master  inertial  navigation  system  (INS)  which  provides  higher 
accuracy  than  the  GGSS'  stabilized  platform  components  can 
provide  (Chapter  5).  A  diagram  illustrating  the  use  of  these 
sensors  as  GGSS  aides  is  presented  in  Fig.  1.2-1. 

Chapter  5  summarizes  only  the  key  master  INS  analysis 
results;  the  complete  technical  details  are  presented  in  Ref.  1 
The  significant  conclusions  of  the  aided-airborne  study  and 
recommendations  of  areas  to  be  pursued  further  are  presented 
in  Chapter  6.  Additional  support  material  and  methodologies 
requiring  extensive  theoretical  development  are  presented  in 
separate  appendices.  Appendix  A  contains  i  GPS  description 
and  discussion  of  anticipated  navigation  accuracies  in  support 
of  operational  surveying  with  the  GGSS;  Appendix  B  presents 
a  derivation  of  the  doppler-shi f t  observable,  an  important 
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measurement  for  moving-receiver  radio  interferometry  to  GPS; 
Appendix  C  discusses  atmospheric  signal  propagation  delay,  a 
significant  error  source  in  applications  of  radio  interferom¬ 
etry;  and  Appendix  D  presents  the  state-space  representation 
of  a  third-order  Markov  process  which  is  used  for  atmospheric 
refractivity  modeling. 


2. 


CONVENTIONAL  GPS-AIDING  AND  BACK-UP 
MODE  ANALYSIS 


The  NAVSTAR  Global  Positioning  System  (GPS)  is  a 
satellite-based,  all-weather,  continuous  radio  navigation 
system  that  provides  accurate  position,  velocity,  and  time 
information  to  users  anywhere  on  or  near  the  earth's  surface. 
When  the  entire  18-satellite  constellation  is  in  place  and 
operational  (late  1988),  conventional  pseudoranging  techniques 
are  more  than  adequate  for  suporting  operational  GGSS  missions. 
As  a  minimum  configuration,  the  current  constellation  of  six 
satellites  will  be  maintained  throughout  the  airborne  GGSS 
test  operations  scheduled  for  the  fall  of  1986.  Despite  the 
limited  viewing  times  associated  with  this  current  constella¬ 
tion,  accurate  position- fixing  and  long-distance  navigation 
using  GPS  have  been  demonstrated  successfully  (Ref.  2).  For 
the  GGSS  application,  satisfactory  GPS-aiding  performance  will 
be  realized  during  testing  through  the  use  of  two  back-up 
modes.  The  first  back-up  mode  involves  incorporating  an  atomic 
clock;  the  second  involves  adding  a  clock  and  aLtimeter.  The 
back-up  modes  are  necessary  to  extend  the  intervals  of  usable 
GPS  coverage  (two  or  more  satellites  in  a  favorable  orientation) 
to  a  sufficient  amount  for  GGSS  test  sorties. 

2.1  SYSTEM  DESCRIPTION  AND  GGSS  MECHANIZATION 

GPS  consists  of  three  major  segments:  1)  a  space 
segment  containing  satellites  which  transmit  radio  signals; 

2)  a  control  segment  made  up  of  ground-based  < quipment  to 
monitor  the  satellites;  and  3)  a  user  segment  comprising 


equipment  which  passively  receives  and  converts  the  satellite- 
transmitted  signals  into  positioning  and  navigation  information. 
The  first  and  third  segments  are  most  relevant  for  this  analy¬ 
sis  and  are  discussed  further  in  Appendix  A.  Additional  details 
are  presented  in  Refs.  3-6. 

2.1.1  GGSS  Mechanization 

As  curren  ly  configured,  the  GGSS  includes  a  TI  4100 
GPS  Navigator,  an  atomic  clock,  and  a  barometric  altimeter 
(see  Fig.  1.1-1).  This  equipment  makes  possible  the  use  of 
the  GPS  back-up  mo  ies  when  only  three  or  two  satellites  are 
visible  during  GGS d  testing  and  initial  survey  operations.  A 
ROLM  MSE-14  milita  ized  computer  controls  the  real-time  data 
flow  and  processing  associated  with  the  GGSS.  A  brief  descrip¬ 
tion  of  the  TI  4100  Navigator  is  presented  below,  followed  by 
a  discussion  of  the  navigation  accuracies  required  for  the 
airborne  GGSS  application. 

The  TI  4100  NAVSTAR  Navigator  consists  of  an  antenna 
assembly,  a  receiver,  and  a  control/display  unit.  'The  antenna 
accommodates  an  omnidirectional  pattern  and  can  receive  both 
the  and  GPS  trequencies  .  The  receiver  is  designed  to 
process  the  signals  from  any  four  visible  GPS  satellies  and 
can  operate  in  a  "degraded"  navigation  and  positioning  mode 
when  a  user  indicates  that  only  three  or  two  satellites  are 
visible .  The  unique  single-channel  hardware  multiplexing 
feature  of  this  receiver  allows  the  and  signals  of  up  to 
four  satellites  to  be  sampled  within  each  20  msec  interval. 

In  a  sampled  d.  ta  sense,  all  the  satellite  signals  are  tracker, 
continuously  as  though  multiple  hardware  receiver  channels  are 
present  and  thus,  interchannel  measurement  bias  is  avoided. 

*GPS  transmits  on  two  carrier  frequencies,  denoted  and  L.;  . 

These  frequencies  are  1575.42  MHz  and  1227.60  MHz,  respectively 


The  output  data  provided  by  the  TI  4100  receiver 
include  the  following : 

•  and  P-code  state  (pseudorange) 

•  and  pseudorange-rate  (P-code) 

•  and  L2  (doppler-induced)  carrier  phase 

(whole  cycles  and  fractions) 

•  Instantaneous  L,  and  Ln  doppler  frequency 

shift  1  1 

•  Averaged  line-of-sight  satellite  acceleration 

•  L-^  C/A-code  state  (pseudorange) 

•  L-^  pseudorange-rate  (C/A-code) 

•  Signal- to-noise  measures  and  other  data 
quality  indicators. 

rhese  data  are  provided  simultaneously  for  four  satellites  at 
three-second  intervals.  Additional  characteristics  and  attri¬ 
butes  of  the  TI  4100  are  presented  in  Ref.  7;  a  comparison  of 
the  characteristics,  capabilities,  and  performance  of  several 
GPS  receivers  currently  available  is  provided  in  Refs.  8  and  9 

2.1.2  GGSS  Navigation  Accuracy  Requirements 


The  core  of  the  GGSS  mechanization  consists  of  three 
gravity  gradiometer  instruments  (GGIs)  mounted  on  a  three-axis 
platform  equipped  with  inertial  sensors  (gyros  and  accelerometers) 
to  provide  local-level  stabilization.  This  mechanization  is 
depicted  in  Fig.  2.1-1.  Extraction  of  precise  gravity  gradients 
requires  1)  accurate  positioning  and  2)  known  orientation  of 
each  GGI .  These  requirements  levy  performance  constraints  on 
the  inertial  platform  to  maintain  misregist ration  and  misresolut ion 
at  acceptable  levels.  Of  particular  interest  in  this  discussion 
is  the  adequacy  of  the  planned  primary  and  >ack-tp  positioning 
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Figure  2.1-1  Schematic  of  Inertially  Stabilized 

Moving-Base  GGSS 

system  in  the  presentee  of  misregistration  due  to  navigational 
position  errors  (i.e.,  measured  gradients  assigned  to  the  wrong 
location  on  the  surface  of  the  earth).  Based  on  previously 
derived  misregistration  sensitivities  for  an  airborne  GGSS 
survey  (Ref.  10),  navigation  accuracies  near  100  m  rms  are 
required  to  avoid  introducing  significant  gradient  misregistra¬ 
tion.  The  discussion  which  follows  quantifies  the  capability 
of  conventional  GPS-aiding  (both  primary  and  back-up  modes)  to 
satisfy  this  navigation  accuracy  requirement. 


2.2  ANALYSIS  AP : ’ROACH  AND  RESULTS 
2.2.1  Approach 

The  analysis  approach  described  herein  is  based  on  a 
simulation  of  the  applicable  GPS  constellation.  Both  the  planned 
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18-satellite  and  the  current  6-satellite  conf igi  rations  are 
considered.  The  technical  details  associated  wi  h  simulating 
a  particular  constellation  are  included  in  Appendix  A.  The 
discussion  in  this  section  emphasizes  the  following  two  points: 
1)  quantifying  the  number  of  GPS  satellite  signals  which  will 
be  available  and  visible  to  the  airborne  GGSS  based  on  the 
orbital  relationships  introduced  in  Appendix  A,  and  2)  assess¬ 
ing  the  associated  measurement  geometry  and  r«sultant  GGSS 
positioning  accuracy. 

Signal  availability  for  radio  communications  systems 
generally  involves  consideration  of  transmitting  power  levels, 
transmitting  antenna  patterns,  the  distance  from  transmitter 
to  receiver,  interference  effects,  receiver  antenna  gains, 
etc.  For  near-earth  GPS  users,  in  particular  for  the  GGSS, 
most  of  these  factors  are  not  applicable  (Ref.  11).  GPS  satel¬ 
lite  orbits  are  so  high  above  the  earth  that  differences  in 
transmitter- to-receiver  distance  are  very  slight  for  near-earth 
users.  Transmitting  antenna  beamwidth  is  adequate  to  cover 
the  entire  near-earth  region  in  the  main  lobe.  GPS  signals, 
having  wavelengths  of  about  19  cm  at  the  L-^  frequency,  travel 
in  a  line-of-sight  fashion  and  are  blocked  by  objects  having 
dimensions  comparable  to  a  wavelength  or  larger  (e.g.,  the 
earth  itself  or  other  obstructions  near  the  receiving  antenna). 
Transmitting  power  levels  and  processing  gains  are  such  that 
C/A-code  acquisition  will  be  possible  for  any  signal  coming 
from  a  GPS  satellite  visible  to  any  part  of  the  receiving 
antenna  pattern  with  a  gain  of  about  -4  dB  or  more.  Thus,  GPS 
signal  availability  for  the  GGSS  is  equivalent  to  satellite 
visibility. 

There  are  four  factors  generally  involved  in  satellite 
visibility: 
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•  The  GPS  satellite  constellation,  i.e., 
the  number  of  operational  in-orbit  satel¬ 
lites  and  their  orbit  relationships 

•  "'he  GGSS  host  aircraft  flight  pattern 
details 

•  The  earth  as  a  signal  block 

•  The  receiving  antenna  coverage. 

The  significance  of  each  of  these  factors  is  discussed  in  Ap¬ 
pendix  A,  along  with  the  navigation  and  positioning  accuracies 
anticipated  from  the  complete  18-satellite  GPS  constellation. 
Key  analysis  results  for  the  six-satellite  constellation  are 
described  in  the  pages  which  follow. 


2.2.2  Six-Satellite  Analysis  Results 

As  pointed  out  in  Appendix  A,  it  is  anticipated  that 
the  GPS  constellation  will  consist  of  only  six  satellites  during 
airborne  GGSS  nesting  scheduled  for  the  fall  of  1986.  This 
condition  has  major  implications  on  the  satellite  visibility 
patterns  which  will  occur  for  the  GGSS  tests  in  North  Texas. 
Figure  2.2-1  illustrates  a  typical  coverage  pattern  which  will 
exist  for  North  Texas  during  a  fall  1986  day.  Since  the  satel¬ 
lites  appear  about  four  minutes  earlier  in  the  sky  each  day, 
the  times  of  adequate  visibility  will  shift  accordingly,  (e.g., 
from  early  morning  to  late  evening  of  the  previous  day,  etc), 
but  the  amount  of  continuous  coverage  will  not  change  appreciably 

Table  2.2-1  quantifies  the  length  of  continuous  cover¬ 
age  associated  with  the  case  shown  in  Fig.  2.2-1.  The  actual 
length  of  usable  coverage  (i.e.,  suitable  for  satisfying  the 
GGSS  navigation  rms  accuracy  goal  of  100  m)  is  determined  by 
examining  the  geometry  associated  with  both  the'  nominal  and 
aiding  solutions  (see  Fig.  2.2-2).  The  aiding  solution  is 
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Figure  2.2-1  Predicted  GPS  Coverage  for  North  Texas  During 

a  Typical  Fall  1986  Day 


TABLE  2.2-1 

SUMMARY  OF  GPS  TYPICAL  COVERAGE  WHICH 
WILL  OCCUR  DURING  GGSS  TESTING* 


MINIMUM  NUMBER 
OF  VISIBLE  SATELLITES 


HOURS  OF 

CONTINUOUS  COVERAGE 
(PRIMARY/SECONDARY  INTERVALS) 


4 . 8/0 . 0 
7 . 0/  0 . 2 
8 . 5/  4 . 5 


rWhen  only  six  satellites  are  available 
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NOMINAL 

SOLUTION 


MEAN 

-  6  3 

STANDARD 

DEVIATION 

=  4.1 

AIOING  ' 
'SOLUTIONS 


3  4  5  6  7  8  9  10  11  12  13  14  15  16  17  18 

TIME  (hr) 


F igure  2.2-2 


PDOP  Characteristics  Associated  with 
Predicted  GPS  Coverage  During 
GGSS  Testing 


accomplished  using  an  atomic  clock  and  barometric  altimeter 
whose  outputs  are  two  and  five  times,  respectively,  less  accur 
ate  than  the  corresponding  quantities  which  would  be  derived 
from  the  nominal  GFS  solution  if  four  satellites  were  visible. 
Aiding  with  a  sufficiently  accurate  atomic  clock  extends  the 
usable  coverage  by  about  55%.  Similarly,  aiding  with  a  suffi¬ 
ciently  accurate  ctock  and  altimeter  provides  an  additional 
24%  of  coverage  during  the  primary  interval  (hours  0  to  12  in 
Fig.  2.2-1)  for  a  total  of  about  8.5  hours.  This  length  of 


-ji  •"  ir»  \ 


-  w  *1  -,» 


time  is  adequate  for  completing  a  daily  50-km  by  50-km  test 
sortie  at  400  km/hr.  In  fact,  four  additional  hours  of  ade¬ 
quate  navigation  are  possible  by  exploiting  the  secondary 
interval  of  GPS  coverage  (hours  12  to  24  in  Fig.  2.2-1).  Note 
that  for  a  typical  flight  pattern,  such  as  that  described  in 
Appendix  A  (Fig.  A.  3-1),  7.5  min  are  needed  to  fly  each  data 
track  and  thus,  2.75  hr  are  required  to  collect  all  the  data. 
The  remainder  of  the  available  time  can  be  expended  in  the 
turns  and  travel  to  and  from  the  calibration  site. 

Also  apparent  in  Fig.  2.2-2  are  two  relatively  short 
(e.g.,  less  than  10  rain)  intervals  when  the  PDOP  exceeds  the 
projected  maximum  value  for  the  airborne  CGSS  application. 
Recall  from  Appendix  A  (Eq.  A. 1-3)  that  the  maximum  PDOP  value 
is  obtained  by  dividing  the  three-dimensional  GGSS  rms  position 
error  goal  (100  m)  by  the  rms  pseudorange  error  (7  m).  Both 
of  these  intervals  of  unfavorable  geometry  are  associated  with 
the  aiding  solutions  and  occur  near  the  end  or  beginning  of 
the  visibility  regions.  Since  these  intervals  are  predictable, 
the  GGSS  testing  schedule  can  be  established  to  appropriately 
avoid  them,  as  necessary.  In  addition,  note  that  the  GGSS 
inertial  platform  is  capable  of  providing  about  5  min  of  ade¬ 
quate  unaided  operation  (Ref.  1). 

2.3  TECHNICAL  RISK  AND  RECOMMENDATIONS 
2.3.1  Technical  Risk 

Conventional  pseudoranging  to  the  NAVSTAR  GPS  will 
provide  the  currently  configured  airborne  GGSS  with  adequate 
navigation  and  positioning  information  during  testing  and 
future  survey  operations.  Nevertheless,  some  technical  and 
program  risk  is  still  associated  with  GPS.  The  key  risk  is 
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that  GPS  is  not  yet  fully  operational.  Fifteen  of  the  planned 
satellites  must  still  be  fabricated,  launched,  and  placed  in 
the  proper  orbit.  Furthermore,  these  satellites  must  remain 
operational  as  long  as  expected  (recall  that  clock  proolems 
contributed  to  the  failure  of  an  early  satellite)  and  the 
position  of  their  orbits  must  be  known  precisely.  Clock  sta¬ 
bility  and  ephemeris  prediction  error  sources  are  the  major 
contributors  to  the  computed  user- to- satellite  rms  ranging 
error  of  seven  meters  rms. 

Another  ar>a  of  risk  is  the  current  status  of  GPS 
user  equipment  in  particular  the  receiver  technology.  Al¬ 
though  the  TI  4100  receiver  has  been  tested  successfully  on 
numerous  occasions  (Refs.  7,  9,  and  12),  this  technology  is 
relatively  new  and  additional  tests,  particularly  in  the  two- 
and  three-satellite  "backup  modes",  are  in  order  (see  below). 

2.3.2  Reconunei  idations 

Based  on  the  analysis  results  presented  herein  and 
published  material,  the  GGSS  Program  can  confidently  proceed 
with  its  planned  use  of  conventional  GPS-aiding.  However, 
to  assure  successfu.  GGSS  testing,  mechanization  of  the  three- 
satellite  and  two-satellite  back-up  modes  is  exceedingly  im¬ 
portant.  In  addition,  to  minimize  the  risk  associated  with 
the  TI  4100  receiver,  extensive  testing  (especially  in  a  dy¬ 
namic  environment  ind  in  tne  degraded  modes)  is  needed  to 
fully  understand  tie  performance  of  the  receiver  hardware 
and  asssociated  soltware.  Finally,  in  anticipation  of  more 
stringent  navigation  and  positioning  requirements  for  the 
GGSS  in  surface  vehicle  applications,  differential  GPS  tech¬ 
niques  should  be  investigated. 
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RADIO  INTERFEROMETRY 


Radio  interferometry  uses  to  advantage  the  wave- like 
nature  of  radio  waves  to  measure  accurately  the  propagation 
path  length  differences  of  two  radio  signals  transmitted  from 
the  same  source.  The  application  of  doppler-shift  and  inter¬ 
fere  metric  phase  analysis  to  GPS  data  collected  in  a  moving- 
receiver  scenario  will  yield  positioning  accuracies  of  a  few 
parts  per  million  (ppm)  for  geodetic  ground  surveys  and  tens  of 
ppm  for  airborne  surveys.  Thus,  radio  interferometry  offers  a 
more  than  adequate  backup  positioning  capability  during  GGSS 
testing  and  survey  operations.  The  potential  for  measuring 
acceleration  to  better  than  one  milligal  accuracy  for  surface 
vehicle  applications  is  also  attractive. 

Data  for  radio  interferometer-based  measurements  are 
obtained  by  employing  the  NAVSTAR  carrier  frequencies  and  thus, 
no  access  to  the  GPS  codes  is  required.  In  addition,  the  GPS 
receiver  included  in  the  nominal  GGSS  configuration  is  uniquely 
qualified  to  provide  the  appropriate  data.  Radio-meteorological 
sensors  and  minor  modifications  to  the  current  GGSS  test  plan 
are  necessary  for  system  implementation.  The  very  accurate 
distance  and  velocity  measurements  possible  using  interferometry 
can  also  be  used  to  augment  a  gravimeter-aided  GGSS  (see  Chap¬ 
ter  4),  thereby  realizing  a  more  accurate  measure  of  gravity. 
The  following  sections  explictly  address  TASC's  analysis  of 
moving-receiver  radio  interferometry. 


3.1  SYSTEM  DESCRIPTION  AND  MECHANIZATION  DETAILS 


Very  high  accuracies  have  been  demonstrated  repeatedly 
during  the  measurement  of  fixed  inter- recei- -er  distances  using 


GPS-based  radio  interferometry.  Phase  observations  recorded 
during  geodetic  surveys  have  resulted  in  accuracies  of  one  ppm 
(or  better)  for  relative  distances  up  to  several  thousand  kilo¬ 
meters  (Refs.  12-15).  However,  the  application  of  radio  inter¬ 
ferometry  to  a  moving- receiver  scenario  is  a  new  area  of  research. 
TASC's  findings  indicate  that  this  new  measurement  technique 
holds  promise  for  the  GGSS  program  as  well  as  the  entire  field 
of  geodesy. 


The  us*  of  radio  interferometry  in  a  dynamic  environment 
is  illustrated  in  Fig.  3.1-1.  At  both  receiver  locations, 
periodically  synchronized  atomic  clocks  are  used  to  time  mark 
the  GPS  signal  phase  and  doppler-shi f t  observations  for  subsequent 
post-mission  data  reduction.  In  addition  to  generating  very 
accurate  position  data  in  support  of  geodetic  survey  missions, 
this  system  car  ser^e  to  supplement  post-mission  reduction  of 


Figure  3.1-1 


Summary  of  Moving-Receiver  Radio 
Interferometry  Concept 


host  vehicle  inertial  navigation  system  (INS)  data  by  providing 
an  accurate  external  velocity  reference.  Furthermore,  by  sub- 
tracting  host  vehicle  INS  and  interferometrically  derived 
accelerations,  one  can  obtain  a  measurement  of  gravity.  This 
is  theoretically  appealing  since  the  accelerometers  in  the  INS 
actually  measure  only  the  difference  between  vehicle  and  gravi¬ 
tational  accelerations  whereas,  using  interferometry,  vehicle 
accelerations  can  be  measured  directly.  Analysis  results, 
presented  in  Section  3.4.2,  quantify  the  accuracy  associated 
with  the  measurement  of  these  navigation  quantities. 

3.1.1  Mechanization  Details 

This  section  addresses  the  minor  modifications  to  the 
GGSS  mechanization  which  are  required  for  radio  interferometry 
to  provide  an  appropriate  backup  positioning  capability  during 
GGSS  test  and  survey  operations.  Depending  on  the  nature  of 
the  gravity  field,  misregistration  sensitivity  studies  indicate 
that  positioning  accuracies  near  10  m  rms  may  be  needed  for 
gravity  surveys  performed  aboard  a  surface  vehicle  (Ref.  10). 
Until  appropriate  testing  procedures  are  carried  out  (e.g., 

Ref.  16)  and  the  collected  data  are  analyzed,  this  stringent 
accuracy  requirement  cannot  be  confirmed.  lowever,  as  noted  in 
Section  2.2  and  Appendix  A  (Section  A. 3),  positioning  accuracies 
near  10  m  rms  will  be  difficult  to  attain  using  standard  GPS 
pseudoranging  techniques.  For  these  reasons,  Bell  Aerospace  and 
7ASC  have  considered  alternative  mechanizations  for  the  surface 
vehicle  implementation.  Note  that  this  issue  is  not  important 
for  airborne  scenarios  in  which  positioning  accuracies  near  100  ra 
rms  are  needed  (Ref.  10).  Accuracy  of  100  m  is  achievable  with 
the  nominal  mechanization  (Sections  2.2  and  A. 3). 


Inertial  Navigation  System. 
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Bell  Aerospace  has  responded  to  the  possible  need  for 
high  positioning  accuracy  during  surface  applications  of  the 
GGSS  by  incorporating  a  "fifth-wheel"  dead-reckoning  mechaniza¬ 
tion  to  the  surface  vehicle  design.  However,  uncertainties 
remain  in  the  fifth-wheel  approach,  particularly  in  regard  to 
the  treatment  of  cross-track  errors. 

As  an  alternative  back-up  approach,  TASC  examined  the 
use  of  GPS-based  radio  interferometry.  This  positioning  system 
requires  modest  augmentation  of  the  current  GGSS/GPS  mechaniza¬ 
tion.  TASC's  analysis  examined  the  ability  of  an  interferometry 
aided  GGSS  to  provide  improved  positioning  accuracies  and  more 
precise  gravity  disturbance  measurements,  compared  with  the 
nominal  GGSS  mechanization  discussed  earlier  in  this  report. 

The  present  GGSS  baseline  configuration  already  include 
a  GPS  receiver,  atonic  clock,  and  data  recording  system.  Thus, 
the  only  new  hardware  needed  at  the  moving  receiver  would  be  a 
radio-meteorological  sensor  to  determine  atmospheric  refractiv- 
ity  variations  on  a  continuous  basis.  In  addition,  an  appropri¬ 
ate  GPS  receiver,  data  recording  system,  refractivity  sensor, 
and  atomic  clock  would  be  required  at  a  fixed-site  location. 

3.1.2  Sensor  Performance 

The  performance  of  any  system  depends  on  the  accuracy 
of  its  consituent  sensors.  For  radio  interferometry,  the  nec¬ 
essary  equipment  consists  of  GPS  receivers  and  atmospheric 
refractivity  monitors.  The  dynamic  environment  further  con¬ 
strains  hardware  to  be  rugged  and  capable  of  maintaining  cal¬ 
ibration  for  the  duration  of  a  survey  sortie  (i.e.  ,  up  to 
8  hours).  These  constraints  limit  the  equipment  available  to 
support  applications  of  moving- receiver  radio  interferometry. 


GPS  Receivers  -  The  receiver  is  required  to  track 
the  carrier  frequency  of  multiple  satellites  in  a  dynamic  en¬ 
vironment;  i.e.,  the  receiver  used  must  maintain  GPS  signal 
tracking  through  doppler  shifts.  As  a  result,  the  hardware 
chosen  must  have  a  wide-band  input  channel.  Neither  the  SERIES- 
ISTAC  receivers,  models  1991  and  2002  (Ref.  17),  nor  the  Macro¬ 
meter  receivers,  models  V-1000  and  II  (Ref.  18)  are  capable  of 
this  since  both  are  designed  as  stationary,  narrow-band  receivers 

The  primary  requirement  of  a  GPS  receiver  for  moving- 
receiver  radio  interferometry  is  to  maintain  sufficient  phase 
accuracy.  Current  technology  dictates  that  this  requirement 
can  only  be  satisfied  by  employing  a  single-channel  multiplexing 
design.  Multichannel  receivers,  such  as  the  SERCEL  TR  55  or 
the  Wild-Magnavox  WM-101,  produce  phase  delay  differences  for 
each  hardware  channel  due  to  the  different  electrical  paths 
encountered  by  each  satellite  signal.  Attempts  have  been  made 
to  compensate  for  these  errors;  however,  these  compensations 
are  still  incapable  of  delivering  the  high  accuracy  required 
for  geodetic  survey  applications  (Ref.  7). 

The  Texas  Instruments  TI  4100  GPS  receiver,  presently 
configured  for  the  baseline  GGSS  system,  is  uniquely  qualified 
to  support  applications  of  moving-receiver  interferometry. 

The  TI  4100  is  the  only  receiver  currently  designed  with  single¬ 
channel  multiplexing.  This  capability  allows  both  the  GPS  L-^ 
'1.5754  GHz)  and  (1.2276  GHz)  frequency  signals  of  up  to 
four  satellites  to  be  sampled  within  each  20-msec  interval 
(Ref.  7).  In  addition,  thj.s  receiver  is  designed  with  a  48-bit 
software  accumulator,  sufficient  to  retain  the  temporal  rela¬ 
tionship  between  carrier  phase  and  time,  as  measured  by  a 
local  clock  for  an  entire  satellite  pass.  The  accumulated 
phase  (whole  cycle  and  fraction)  is  provided  to  the  user  at 
’.0-  and  1.2-second  intervals  via  a  9600-baud  RS-232  port.  The 


phase  is  derived  with  respect  to  either  the  internal  receiver 
oscillator  (with  20-nsec  accuracy)  or  an  external  user- suppl ied 
clock  which  can  augment  the  receiver.  (The  time  can  also  be 
recorded  along  with  the  signal  phase  for  post-mission  process¬ 
ing.  )  In  one  demonstration  conducted  by  the  Geodetic  Survey 
Division,  Surveys  and  Mapping  Branch  of  Energy,  Mines,  and 
Resources  of  Canada,  recorded  phase  observations  from  two 
TI  4100s  were  used  to  recover  inter-receiver  distances  up  to 
1C2  km  in  length  with  accuracies  of  one  ppm  in  all  three  co¬ 
ordinates  (Ref.  19). 


In  addition,  the  TI  4100  outputs  the  "instantaneous" 
doppler  frequency  shift  observed  for  each  satellite  signal. 
This  measurement  is  the  result  of  endpoint  averaging  two  fre¬ 
quency  shifts  160  msec  apart  for  both  the  and  I^  signals. 
Recording  of  both  the  phase  and  doppler- shi ft  outputs  from  the 
TI  4100  can  be  readily  accomplished  using  the  planned  16  or 
0.5  Hz  sampling  rate  of  the  GGSS  data  acquisition  and  monitor¬ 
ing  system.  A  summary  of  the  attributes  and  limitations  of  the 
receivers  reviewed  for  this  study  is  presented  in  Table  3.1-1. 

Atmospheric  Monitors  -  The  largest  source  of  error 
for  interferoraetrically  derived  quantites  is  due  to  signal 
delay  in  the  troposphere.  In  order  to  measure  inter-receiver 
distances  using  radio  interferometry,  the  time  delay  between 
the  receipt  of  the  same  signal  at  two  separate  receivers  must 
be  recorded.  This  time  delay  (which  causes  the  receivers  to 
be  "out  of  phase"  with  each  other)  is  affected  by  the  differing 
refractive  properties  of  the  atmosphere  along  the  radio  path 
at  the  two  receiver  locations. 


Atmospheric  refractivity  is  influenced  by  both  elec¬ 
tromagnetic  conditions  in  the  ionosphere  and  meteorological 
conditions  in  the  troposphere.  (These  affects  are  described 
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TABLE  3.1-1 

CAPABILITY  OF  AVAILABLE  GPS  RECEIVERS 
TO  SUPPORT  MOVING-RECEIVER  RADIO  INTERFEROMETRY 


RECEIVER 

ATTRIBUTES 

LIMITATIONS 

ISTAC-SERIES 
(1991  &  2002) 

■■ 

Codeless,  lightweight 

Narrow  bandwidth,  stationary 
applications 

Macrometer 

(V-1000  &  II) 

Six-satellite  tracking, 
codeless 

Narrow  bandwidth,  stationary 
applications 

SERCEL  TR  55  !  J 

(  Accurate  dynamic  position- 
Wild-Magnavox  (  ing  via  codes 

WM- 101  !  ' 

_ _ _ J _ 

Multiple  hardware  channels, 
limited  test  results 

Texas  Instruments 
TI  4100 

Single-channel  multiplexing, 
dynamic  environment  design, 
two-  and  three-satellite 
degraded  mode  design 

Limited  test  results 

in  Appendix  C.)  Refractivity  is  measured  in  N-units  which 
are  defined  as 

N  =  (c/v  -  1)  x  106  (3.1-1) 

Q 

where  c,  the  vacuum  speed  of  light,  is  2.9978  x  10  m/sec 
and  v  is  the  atmospheric  (refracted)  speed  of  light.  Thus, 
a  change  of  one  N-unit  corresponds  to  a  nanosecond  delay  for 
typical  atmospheric  signal  paths  (i.e.  about  60  km  of  tropo¬ 
sphere).  Measurement  of,  and  correction  for,  atmospheric 
oropagation  delays  are  described  in  Appendix  C. 

Ionospheric  delays  can  be  compensated  by  modeling 
the  upper  atmosphere  or  by  using  dual- frequency  radio  receivers. 
The  former  approach  reduces  signal  delays  by  50  to  70  percent 
on  an  rms  basis)  compared  with  uncorrected  leasurements ; 

Jual - frequency  receivers  reduce  signal  delays  by  better  than 
-•0  percent  on  an  rms  basis. 


I 
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Tropospheric  refractivity  can  be  measured  using  three 
different  atmospheric  sensors  (i.e.,  radiosondes,  refractom- 
eters,  and  radiometers).  Radiosondes ,  generally  employed  for 
weather  forecasting,  have  limited  sensitivity  and  slow  response 
to  the  varying  atmospheric  conditions  expected  during  applica¬ 
tions  of  moving- receiver  interferometry.  Microwave  ref ractom- 
eters  measure  changes  in  cavity  resonance  and  are  most  frequently 
used  to  record  refractivity.  Radiometers  measure  thermal 
atmospheric  emission  and  have  only  recently  been  developed  and 
tested  for  recording  line-of-sight  propagation  delay  times. 
Although  still  in  an  experimental  stage,  portable  radiometers 
are  expected  soon. 


The  three  radio-meteorological  sensors  examined  for 
use  in  moving-receiver  radio  interferometry  are  listed  in 
Table  3.1-2.  Subsequent  analysis  results  which  are  reported 
herein  reflect  the  availability  of  ref ractometers  to  monitor 
atmospheric  variations. 


TABLE  3.1-2 

EVALUATION  OF  RADIO-METEROLOGICAL  SENSORS 


SENSOR  FAMILY 

PROPERTY  MEASURED 

LIMITATION 

RMS 

ACCURACY 

(N-UNITS) 

Weather  radiosonde 

Air  pressure,  humidity, 
temperature 

Slow  response  to 

refractivity 

variations 

1.0 

Ref ractometer 

Refractivity  variations 

Initial  calibration 
required 

0.01 

Dual- frequency 
radiometer 

H^O  thermal  atmospheric 

emission  at  22.2  and 

31.4  GHz 

Still  in  experimental 
stage,  not  portable 

0.0001 

3.2 


DATA- REDUCT ION  APPROACH 


3.2.1  Interferometric  Analysis 

Radio  interferometry  requires  recordin ?  of  the  time 
delay  between  the  receipt  of  the  same  signaL  at  two  separate 
receivers  (see  Fig.  3.2-1).  The  time  delay  determines  the 
amount  of  phase  retardation  present  between  receivers;  i.e., 
at  a  given  instant  in  time,  how  many  carrier  cy:les  exist  be¬ 
tween  the  two  receivers.  This  observed  phase  difference,  A<])  , 
can  be  related  to  the  perpendicular  distance,  d,  from  the 
moving  receiver  to  the  signal  wavefront  plane  through  the 
equation 


d  =  cA0/(nf)  (3.2-1) 


where  f  is  the  GPS  L-^  signal  carrier  frequency  (1.5754  GHz) 
and  n  is  the  index  of  refraction.  A  subsequent  measurement  of 
the  phase  difference  reveals  a  new  value  for  d  and  yields  the 
GPS  signal  plane's  angle  of  incidence,  0,  with  respect  to  the 

a-jioto 


V  -  VEHICLE  SPEED 

At  -  TIME  SAMPLING  INTERVAL 

d  -  PHASE  DIFFERENCE  PERPENDICULAR  DISTANCE 

6  -  SATELLITE  SIGNAL  INCIDENCE  ANGLE 

x  -  HORIZONTAL  INTER-RECEIVER  DISTANCE 

h  -  ALTITUDE 


Figure  3.2-1 


i 


Overview  of  Interferometric  Phase -Di f ference 
Geometry  in  Two  Dimensions 
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vehicle's  direction  of  motion.  The  motion  is  expressed  as  a 
function  of  vehicle  velocity,  v,  and  time  between  phase  measure 
ments,  t.  By  employing  three  or  more  satellites  from  the  GPS 
constellation,  the  inter- receiver  distance  can  be  uniquely 
determined  (within  a  constant  ambiguity  resolved  through  proper 
initialization) . 


The  horizontal  inter- receiver  distance,  x,  can  be 
related  to  the  phase  difference  through  d  and  the  equation 


x  =  (h/tar0)  +  (d/sin0) 
which  is  the  equivalent  to 

x  =  (nfhccs0  +  cA<J>  )/(nfsin0  )  (3.2-2) 

where  h  is  the  survey  vehicle  altitude.  In  terms  of  A$  , 
Eq.  3.2-2  becomes 

A<}»  =  ^  (xsin©  -  hcos0 )  (3.2-3) 


3.2.2  Doppler  Analysis 

An  improved  measurement  of  the  vehicle  velocity  can 
also  be  obtained  from  the  GPS  signals  using  a  doppler- shi f t 
analysis.  As  a  receiver  moves  with  respect  to  an  emitting 
source,  the  received  frequency  is  doppler- shi f ted .  For  radio 
interferometry  to  GPS,  satellite  motion  introduces  two  doppler- 
shift  frequencies  (or.e  at  each  receiver).  These  frequencies 
are  inherently  different  due  to  the  additional  velocity  associ¬ 
ated  with  the  moving  receiver,  and  can  be  related  by  the  equation 


r  1  *  nvcos0/c  i 
1  -  nvcosQ/c  j 


(3.2-4) 
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where 


f  =  the  doppler-shifted  frequency  as  re¬ 
ceived  at  the  fixed-receiver 
f'  =  the  doppler-shifted  frequency  as  received 
at  the  moving  receiver  and  due  to  the 
vehicle  motion. 

The  difference,  Af,  between  the  moving-receiver  frequency 
and  the  fixed-receiver  frequency  represents  the  doppler-shift 
due  to  moving-receiver  velocity  alone  and  car  be  derived  from 
Eq.  3.2-4  as 

af  =  fvngose  j  +  vnegse  1  (3.2-5) 

Details  of  the  square-root  expansion  used  to  obtain  Eq.  3.2-5 
are  presented  in  Appendix  B. 


3.3  APPLICABLE  ERROR  SOURCES  AND  MODELS 

Interferometrically-derived  distance  measurements 
require  knowledge  of  the  difference  between  signal  receiver 
paths.  Hence,  many  of  the  error  sources  common  to  pseudorang¬ 
ing  (single-receiver)  techniques  do  not  affect  dual-receiver 
measurement  integrity.  Error  sources  inherent  to  fixed-site 
radio  interferometry  and  which  are  also  applicable  in  dynamic 
aDplications ,  include  receiver  noise,  emitted  frequency  insta¬ 
bility,  satellite  ephemeris  error,  and  clock  jitter  as  well 
as  previously  discussed  (Section  2.2)  satellite  geometry  con¬ 
straints  (i.e.,  GDOP  effects).  In  addition,  there  are  mechani¬ 
cal  and  environmental  error  sources  which  are  unique  to  a 

*This  doppler  shifting  is  due  to  satellite  motion  relative  to 
the  receivers  and  is  approximately  the  same  for  both  receivers. 


moving-receiver  mechanization.  These  include  relative  timing 
drift,  atmospheric  propagation  errors,  and  multipath  effects. 
The  applicable  error  sources  are  addressed  below  by  category 
(i.e.,  environmental,  mechanical,  measurement). 

3.3.1  Environmental 


Tropospheric  Propagation  Delay  (Horizontal  Variations) 
As  described  in  Appendix  C  (Section  C.2),  effects  in  the  lower 
atmosphere  are  an  important  contributor  to  interferometric 


measurement  uncertainty.  Tropospheric  refractivity  varies 
both  temporally  and  spatially  from  survey  point  to  point. 
Furthermore,  correct,  modeling  of  these  fluctuations  requires 
regional  (horizontal)  atmospheric  variabilities  to  be  accounted 
for.  One  proposed  approach,  known  as  site  diversity  (Ref.  20), 
considers  an  exponential  correlation  between  atmospheric  condi¬ 
tions  at  the  two  receiver  sites.  This  approach  uses  only  one 
radio-meteorological  sensor  and  neglects  refractivity  effects 
due  to  pressure  fronts  and  clouds.  Based  on  site  diversity,  a 
statistical  model  can  be  implemented  using  a  first-order  Markov 
process  to  characterize  the  correlation  effect.  To  study  inter¬ 
ferometry  measurement  accuracies,  empirical  data  were  used  to 
determine  appropriate  parameters  for  the  Markov  process.  A 
correlation  distance  of  30  km  and  an  rms  error  of  3  N-units 


were  selected. 


An  alternative  approach  to  site  diversity  requires 
the  use  of  radio-meteorological  sensors  at  both  receiver  loca¬ 
tions.  Errors  associated  with  this  approach  are  a  function  of 
the  particular  sensor  used  (see  Table  3.1-2)  and  can  be  statis¬ 
tically  modeled  using  a  white-noise  process.  Appropriate 
parameters  for  this  model  are  an  rms  error  of  0.01  N-units  and 
an  averaging  time  of  0.1  sec. 


3-12 


Tropospheric  Propagation  Delay  (Vertical  Variations)  - 
Since  vertical  variations  (fluctuations)  in  tropospheric  refrac- 
tivity  increase  with  signal  path  length  (see  Appendix  C,  Sec¬ 
tion  0.2),  use  of  a  spherical  layer  representation  of  the 
atmosphere  above  the  earth  offers  an  advantageous  modeling 
geometry.  In  each  layer,  the  refractivity  fluctuations  vary 
as  the  cosecant  of  the  elevation  angle  0.  ’hese  variations 
can  be  statistically  modeled  by  a  white-noise  process  with  an 
rms  value  of  0.4  arc  min  (30-deg  user-to-satellite  elevation 
angle)  for  six-second  averaging  times  (Ref.  20). 

Multipath  Effects  -  Moving-receiver  radio  interferom¬ 
etry  data  can  be  corrupted  by  two  sources  of  multipath  errors. 
The  first  is  due  to  vehicle-induced  scattering  of  the  GPS 
signal,  i.e.,  temporary  signal  reflections  caused  by  adjacent 
structures  or  emanating  from  vehicle  surfaces.  These  reflec¬ 
tions  can  be  statistically  modeled  by  a  random  process  and 
require  further  investigation.  For  analysis  purposes,  this 
error  is  modeled  as  white  noise  with  an  rms  error  of  1  x  10  ^  m 
for  a  2.5-sec  averaging  time.  The  second  source  of  multipath 
propagation  is  a  result  of  atmospheric  ducting  (signal  bounces 
from  ground  to  atmosphere).  Such  reflections  occur  at  low 
satellite  elevation  angles.  These  errors  are  avoided  by  requir¬ 
ing  satellite  elevation  angles  to  be  greater  than  30  degrees. 

Ionospheric  Propagation  Delay  -  The  propagation  delay 
caused  by  the  passage  of  electromagnetic  waves  through  the 
ionosphere  is  discussed  in  Appendix  C  (Section  C.3).  For 
interferometric  applications,  this  signal  delay  contributes 
measurement  uncertainties  whenever  the  inter-receiver  distances 
are  such  that  the  GPS  signal  travels  along  two  significantly 
different  ionospheric  paths  (i.e.,  distances  greater  than 
1000  km).  Furthermore,  this  delay  can  be  adequately  compen- 
cated  through  the  use  of  dual- frequency  receivers  or  by  appro¬ 
priately  modeling  the  upper  atmosphere  during  data  preprocessing 
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(see  Appendix  C) .  For  this  reason,  effective  compensation  of 
ionospheric  sigral  delay  is  assumed  in  TASC ' s  analysis. 

3.3.2  Mechanical 

Clock  Errors  -  To  minimize  timing  errors  when  record¬ 
ing  signal-phase  differences,  a  moving-receiver  interferometry- 
based  system  requires  that  the  atomic  clocks  at  the  receiver 
sites  be  periodically  synchronized.  For  cesium  time  standards, 
synchronization  can  be  accomplished  by  phase-matching  and  time¬ 
averaging  clock  sire  waves,  resulting  in  correlation  between 
the  time  pieces  of  one  part  in  10^  for  100-second  averaging 
times  (Ref.  21)  After  synchronization,  the  clocks  will  decor¬ 
relate.  The  short-term  loss  of  synchronization  ("uncompensated 
relative  drift" )  is  statistically  modeled  using  a  random-ramp 
analysis.  A  random  ramp  is  typically  used  to  describe  errors 
which  exhibit  time-growing  behavior.  The  growth  rate  of  the 
random  ramp  is  a  random  quantity  with  a  given  probability 
density  (Ref .  22 ) . 

In  addition,  the  clocks  have  individual  "jitter." 
This  error  can  be  characterized  using  a  random  walk  analysis. 

A  random  walk  process  results  when  uncorrelated  signals  are 
integrated.  This  situation  occurs  when  determining  the  cumu¬ 
lative  effect  of  clock  noise.  For  cesium  clocks  (e.g.,  the 
HP  5061A)  ,  the  rms  error  is  1  x  10”^  sec  (Ref.  23);  thus,  in 
six  hours  (a  typical  GGSS  survey  sortie)  the  system  timing 
noise  is  two  nanoseconds.  This  introduces  a  10-cm  error  in 
inter-receiver  distance  measurements. 


Emi tte d  Frequency  Instability  -  Analysis  of  the  oscil¬ 


lators  aboard  the  GPS  satellites  (Ref.  24)  has  demonstrated  a 
stability  of  the  emitted  frequency  to  one  part  in  10^  for 
averaging  times  of  10^  sec.  This  instability  yields  distance 


T5V7r 


measurement  errors  on  the  order  of  10  ra.  Thus,  this  error 
source  can  be  neglected  for  the  accuracy  level  anticipated  for 
the  GGSS  application. 

Altitude  Perturbations  -  Any  variation  in  the  altitude 
of  the  moving  receiver  between  two  sampling  times  will  perturb 
the  observed  signal  incidence  angle,  0.  This  phenomenon  can 
be  seen  in  Fig.  3.3-1.  The  observed  distance,  ASq,  from  the 
satellite  signal  wavefront  to  the  second  measurement  point  B 
can  be  written 


AS  _  =  AS  +  AS„ 
o  t  e 


.  (3.3-1) 


where  the  subscripts  t  and  e  refer  to  the  true  and  error  por¬ 
tions  of  the  observed  distance.  To  a  close  approximation,  the 
error  contribution  to  the  observed  perpendicular  distance  is 
given  by 


-  VEHICLE  SPEED 

-  TIME  SAMPLING  INTERVAL 

-  OBSERVED  PERPENDICULAR  DISTANCE 

-  TRUE  PERPENDICULAR  DISTANCE 

-  ERROR  PERPENDICULAR  OISTANCE 

-  ALTITUDE  perturbation 

-  SATELLITE  SIGNAL  INCIDENCE  ANGLE 


: 

rv\ 


B  \  \AS, 


Figure  3.3-1 


Effect  of  Altitude  Perturbation  on 
Signal  Incidence  Angle 
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.  -  A  A 


(3.3-2) 


AS£  =  <5h/cos6 

Employing  small-angle  approximations,  the  error  for  the  incident 
angle,  60,  can  be  written 

60  =  AS£/( vtcos0 ) 

or  equivalently, 

60  =  6h/(vtcos^0)  (3.3-3) 

The  altitude  perturbation,  6h,  is  a  vehicle-dependent  error 
source  representing  the  amount  a  vehicle  will  "bounce"  away 
from  the  nominal  trajectory.  This  error  source  can  be  statis¬ 
tically  modeled  using  a  white-noise  process  with  an  rms  error 
of  1.0  m  over  a  2.5-sec  averaging  time. 

2 

The  theoretical  effect  of  the  (1/cos  0)  factor  in 
Eq.  3.3-4  is  to  create  a  singularity  at  0  =  90°;  i.e.,  the 
situation  when  a  satellite  is  directly  above  the  two  receivers. 
In  practice,  as  a  satellite  climbs  to  the  overhead  position, 
the  GDOP  degrades  until  another  satellite  is  selected  (yielding 
a  more  favorable  geometry). 

Satellite  Ephemeris  Error  -  Radio  interferometry  to 
GPS  satellites  is  a  generalization  of  quasar-based  stellar 
interferometry.  In  actuality,  the  GPS  satellites  are  not  very 
distant,  compared  with  quasars,  so  that  the  satellite  signals 
arrive  at  the  earth  as  spherical  wavefronts.  To  apply  the 
experience  gained  from  very  long  baseline  interferometry  (VLBI) 
geodesy  to  the  case  of  GPS  signal  sources,  the  actual  spherical 
wavefronts  must  be  reformulated  into  a  plane  wave  with  a  cor¬ 
rection  term  requiring  explicit  range  information.  The  depar¬ 
ture  of  the  GPS  wavefronts  from  plane  waves  is  functionally 
represented  by 


AS  a  b  /r 


(3.3-4) 


where  b  is  the  inter-receiver  distance  and  r  is  the  satellite 
orbit  radius  (r  ~  26  x  10^m) .  The  ephemeris  error  is  propor¬ 
tional  to  As/r  and  leads  to  an  error  of  one  millimeter  for 
inter-receiver  distances  of  1000  km.  This  contributor  to 
navigation  error  can  be  statistically  modeled  by  a  random 
constant  for  inter-receiver  distances  greater  than  1000  km. 

For  scenarios  involving  distances  shorter  than  this,  satellite 
ephemeris  error  can  be  neglected. 

3.3.3  Measurement 

Receiver  Limitations  -  Accurate  measurement  of  the 
phase  difference  and  doppler  shift  between  two  receiver  sites 
ultimately  depends  on  the  GPS  receiver  design.  The  TI  4100 
receiver  incorporates  a  time-sharing  digital  oscillator  for 
multiple  satellite  signal  tracking  and  requires  precise  car¬ 
rier  phase  initialization  at  the  start  of  each  satellite  signal 
sampling  period.  This  carrier  phase  preset  introduces  a  phase¬ 
tracking  rms  error  (Ref.  7)  of  1.525  *  10  ^  cycles  per  cycle 
(averaging  time  of  2.5  sec)  and  is  included  as  measurement 
noise  (Ref.  22)  in  the  interferometry  analysis. 

The  GPS  receiver  is  designed  with  phase  lock  loop 
(PLL)  circuitry.  By  necessity,  the  PLL  tracking  bandwidth 
must  be  wide  enough  to  accept  frequency  shifts  caused  by  the 
relative  motions  of  the  satellite  and  receiver.  The  design  of 
the  receiver  introduces  a  doppler-shif t  measurement  error  of 
0.7  Hz  rms  for  an  averaging  time  of  2.5  sec.  During  the  navi¬ 
gation  analysis,  this  effect  is  modeled  as  measurement  noise. 

Summary  of  Error  Sources  -  A  compilation  of  the  error 
sources  associated  with  moving-receiver  radio  interferometry 
is  presented  in  Table  3.3-1.  For  each  source,  the  error  model 
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TABLE  3.3-1 


COMPILATION  OF  MOVING-RECEIVER 
RADIO  INTERFEROMETRY  ERROR  SOURCES 


7-  '  —o'* 


DESCRIPTION 

MODEL  STRUCTURE/ 
PARAMETERS 

VALUE 

UNITS 

ENVIRONMENTAL 

Tropospheric  Propagation  Delay: 

Horizontal  Variations:  # 

Single-site  monitoring 

Markov  process/ 

Correlation  distance 

30 

km 

RMS  error 

3 

N-uni ts 

Dual-site  monitorin 

White  noise/ 

RMS  error 

0.01 

N -uni t s 

Averaging  time 

0.1 

sec 

Vertical  Variations 

White  noise/ 

RMS  error 

0  .  U 

arc  min 

Averaging  time 

b.O 

sec 

Multipath  Effects: 

Vehic le- induced 

White  noise/ 

o 

o 

V 

-~4 

RMS  error 

m 

Averaging  time 

2.5 

sec 

Ground- induced 

Masked  out  by  elevation 
angle  requirements 

--- 

... 

Ionospheric  Propagation  ielay 

Not  modeled^ 

... 

... 

MECHANICAL 

Clock  Errors: 

Uncompensated  relative  drift 

Random  ramp/ 

1  *  10‘U 

RMS  ramp  slope 

sec./sec 

Compensated  relative  drift 

( c loc  k  jitter; 

Random  walk/ 

RMS  error 

1  *  10‘13 

sec 

Emitted  Frequency  Instability 

White  noise/ 

1  1  10'10 

RMS  error 

Hz 

Averaging  time 

1  *  105 

sec 

Altitude  Perturbation 

White  noise/ 

RMS  error 

1.0 

m 

Averaging  time 

2.5 

sec 

Satellite  Ephemeris  Error 

Not  Modeied1^ 

... 

... 

MEASUREMENT 

Phase-Tracking  Error 

White  noise/ 

RMS  error 

1.525  *  10  3 

eye les/eve  Le 

Averaging  time 

2.5 

sec 

Doppler-shi ft  Measurement  Error 

White  noise/ 

RMS  error 

0.7 

Hz 

Averaging  time 

2  .  5 

sec 

-Not  part  of  nominal  analvsis;  included  for  separate  sensitivity  analyses. 
f Not  applicable  for  inte r - rece 1 ve r  distances  less  than  1000  km. 
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and  associated  parameter  values  which  were  used  in  TASC ' s 
analysis  are  provided  in  the  table.  Discussion  of  study  re¬ 
sults  follows. 

3.4  ANALYSIS  RESULTS 

3.4.1  Simulation  Details 

The  position,  velocity,  and  acceleration  measurement 
capability  of  moving-receiver  radio  interferouetry  was  assessed 
using  covariance  analysis  techniques  and  the  error  models  des¬ 
cribed  in  the  previous  section.  In  this  analysis,  a  coplanar 
base-station,  host  vehicle  trajectory,  and  single-satellite 
flight  path  were  simulated.  Optimal  Kalman  smoothing  tech¬ 
niques  (Ref.  22)  were  applied  to:  1)  determine  the  interfer¬ 
ometer  system's  ability  to  measure  navigatio i  quantities  and 
2)  study  the  sensitivity  of  measurement  accuracy  to  vehicle 
velocity,  atmospheric  refractivity  sensor  availability,  timing 
drift  compensation,  and  vehicle- to-satellite  elevation  angle. 
Both  surface  and  airborne  vehicle  surveying  applications  were 
simulated . 

The  error  dynamics  of  the  moving-receiver  system  are 
represented  by  the  linear  first-order  vector-m;  trix  differen¬ 
tial  equation 

x(t)  =  F(t)  x(t)  +  w(t)  (3.4-1) 

where  x(t)  is  the  system  state  vector,  w(t)  is  the  random 
forcing  function,  and  F(t)  is  the  system  dynamics  matrix. 
Using  this  form,  the  following  expressions  can  be  written  for 
the  vehicle  horizontal  position  (5p),  velocity  (6v) ,  and  accel¬ 
eration  (5a)  errors: 
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where 


6p  = 

6v  = 
6a  = 


-v/n  6n  +  vtanQ  60  +  £v/f  6f 


-v/(nt)  6n  +  vtan0/t  60  +  £v/(ft)  6f  -  v/l  6t 


-v/(nt^)  5a  +  vtan0/t^  60  +  £v/(ft^)  6f  -  v/t^  6 


(3.4-2) 


v  = 


nominal  vehicle  speed  (50  km/hr  -  surface 
van;  450  km/hr  -  aircraft) 


n 

0 

f 

t 

4 

6n 
60 
6  f 
6  t 


nominal  index  of  refraction  (1.00035) 
receiver  to  satellite  elevation  angle  (30  deg) 
nominal  GPS  transmitting  frequency  (1.5754  GHz) 


=  averaged  measurement  interval  (1.2  sec) 


2  2 
c  -  (nvcosQ)  ]/(cncos0) 


tropospheric  refractivity  error 
signal  elevation  angular  error 
emitted  frequency  instability 
clock  errors. 


The  velocity-dependent  term,  4,  and  the  differential  equation 
for  6p  are  derived  in  Appendix  B. 


The  Simula  ion  scenario  incorporates  a  single  GPS 
satellite,  a  TI  4100  receiver  at  a  fixed-site,  and  a  TI  4100 
receiver  aboard  a  moving  vehicle  (surface  van  and  aircraft). 
The  moving  receiver  travels  at  a  constant  velocity  and  alti¬ 
tude  away  from  the  fixed  receiver.  In  addition,  the  flight 
paths  of  the  satellite  and  the  vehicle  are  coplanar  and  inter¬ 
receiver  distances  Less  than  1000  km  were  maintained.  Note 
that  for  distances  larger  than  this,  the  flat-earth  approxima¬ 
tions  implicit  in  t  he  dynamics  equations  (e.g.,  identical 


satellite  frequencies  and  elevation  angles  at  both  receivers) 
are  no  longer  valid.  In  all  simulations,  the  atomic  clocks 
are  considered  synchronized  within  nominal  phase- tracking 
error  limits  at  survey  initialization  (see  clock  error  models 
listed  in  Table  3.3-1).  A  summary  of  the  nominal  scenario 
simulated  in  TASC's  analysis  of  moving- rece Lver  radio  inter¬ 
ferometry  is  presented  in  Table  3.4-1. 

3.4.2  Results 

This  section  describes  the  results  associated  with 
TASC's  study  of  moving-receiver  radio  interferometry  to  support 
moving-base  gravity  field  survey  activities.  In  particular, 
the  relative  distance,  velocity,  and  acceleration  measurement 
capability  of  the  system  for  both  surface  vel icle  and  airborne 
applications  is  addressed.  For  these  analyses,  the  simulation 
explicitly  accounts  for  the  use  of  dual-site  refractivity  mon¬ 
itors  and  presumes  that  post-mission  compensation  for  relative 


TABLE  3.4-1 

NOMINAL  SCENARIO  DETAILS 


ITEM 

DESCRIPTION 

Ha  rdwa re 

Gravity  Sensor 

Clock 

GPS  Receiver 
Radio-Meteorological 
Sensor 

GGSS 

Atomic  (at  fixed  site  and  aboarl  moving  vehicle; 

relative  timing  drift  compensated  for) 

TI  4100  (at  fixed  site  and  aboard  moving  vehicle) 
Microwave  ref ractometer  (at  both  receiver  locations) 

Operating  Conditions 

Vehicle  Type 

Van  Speed 

Aircraft  Speed 

Aircraft  Altitude 
Vehicle  Flight  Path 
Signal  Elevation  Angle 

Surface  van,  aircraft 

50  km/hr 

450  km/hr 

600  m 

Coplanar  with  satellite;  a /ay  from  fixed  receiver 

At  least  30  deg 
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clock  drift  can  be  performed  to  accuracies  already  demonstrated 
by  other  invest igat ions  (Ref.  15).  In  addition,  the  sensitiv¬ 
ity  of  the  mea  jure  aent  accuracy  to  key  operational  elements 
(vehicle  speed,  atmospheric  refractivity  sensor  availability, 
timing  drift  compe lsation ,  and  signal  elevation  angle)  is 
quantified . 

The  navigation  accuracies  associated  with  surface 
vehicle  and  airborre  applications  of  radio  interferometry  are 
presented  in  T  lble  3.4-2.  These  results  indicate  that  for  a 
50-km  survey  track,  position  accuracy  to  better  than  one-meter 
ca:i  be  anticipated;  similarly,  for  a  500-km  track,  less  than 
10  m  error  is  expected.  Results  associated  with  varying  the 
vehicle  speed  up  to  450  km/hr  are  presented  in  Fig.  3.4-1. 
Note  that  for  near  stationary  receivers,  the  positioning  ac¬ 
curacy  approaches  <■  bout  1-2  ppm,  which  is  the  accuracy  level 
typically  reported  from  fixed-site  investigations  (Refs.  13-15) 
These  encouraging  positioning  accuracies  demonstrate  the  ade¬ 
quacy  of  usinp,  this  approach  as  a  back-up  during  GGSS  testing 
particularly  in  surface  vehicle  testing  scenarios  where  strin 
gent  Dositioning  a  :curacies  may  be  required. 


TABLE  3.4-2 

MEASUREMENT  CAPABILITY  OF  RADIO  INTERFEROMETRY 


NAVIGATIC  N 
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0.002 
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mgal 
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Figure  3.4-1  Sensitivity  of  Positior  Measurement 

Accuracy  to  Vehicle  Speed 

Vehicle  speed  has  a  significant  efftct  on  tie  acceler¬ 
ation  measurement  capability  of  radio  inter,  eroniet.y.  This 
effect  is  illustrated  in  Fig.  3.4-2.  For  slowly  moving  surface 
vehicles,  a  theoretical  accuracy  level  on  the  order  of  that 
anticipated  from  the  GGSS  is  anticipated.  The  degradation  of 
acceleration  measurement  accuracy  to  4.4  mgal  f‘>r  the  airborne 
application  is  necessarily  due  to  the  increased  vertical  posi¬ 
tion  and  velocity  uncertainty  (typically,  an  order  of  magnitude 
difference)  at  the  higher  speeds. 

Sensitivity  of  Measurement  Accuracy  to  Atmospheric 
Refractivity  Sensor  Availability  -  The  effect  o'  tropospheric 
refractivity  on  interferometric  measurement  acc  iracy  has  been 
the  subject  of  extensive  error  modeling  (Res.  20  and  25)  and 
hardware  development  (Refs.  26  and  27).  The  signal  delay 
■  aused  by  the  lower  atmosphere  requires  at  least  one  monitor 


' 
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Figure  3.4-2  Sensitivity  of  Acceleration  Measurement 

Accuracy  to  Vehicle  Speed 

(situated  at  either  receiver  location)  to  record  local  atmos¬ 
pheric  variations.  However,  since  a  moving- receiver  applica¬ 
tion  of  radio  interferometry  increases  the  probability  of 
encountering  tropospheric  conditions  significantly  different 
from  fixed-site  applications  (hence,  different  refractive 
indices),  atmospheric  monitoring  at  both  receiver  sites  is 
appropriate.  The  implications  of  single  vs  dual-site  monitor¬ 
ing  on  measurement  accuracy  are  discussed  in  the  paragraphs 
which  follow. 

As  not<d  in  Section  3.1.1  and  summarized  in  Table  3.1-1, 
error  models  for  tropospheric  signal  delay  are  included  in  the 
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analysis.  Two  scenarios  were  simulated  as  part  of  this  sensitiv¬ 
ity  study.  The  first  considered  a  single  microwave  refractometer 
monitoring  local  atmospheric  fluctuations  at  the  fixed-site 
receiver.  A  third-order  Markov  process  (see  Table  3.3-1)  was 
used  to  provide  the  necessary  correlation  between  the  atmos¬ 
pheric  conditions  at  the  two  receiver  sites.  (The  derivation 
of  the  matrices  associated  with  the  state-space  form  of  a 
third-order  Markov  process  is  included  in  Appendix  D.)  The 
second  scenario  simulated  the  case  involving  microwave  refrac- 
tometers  located  at  both  receiver  locations.  (Note  that  this 
is  the  nominal  system  configuration  considered  in  TASC ' s  analy¬ 
ses  as  indicated  in  Table  3.4-1.)  Refractivity  measurement 
inaccuracy  is  modeled  as  white  noise  (see  Table  3.3-1)  in  this 
case.  Other  details  of  the  simulated  scenario  were  unchanged 
from  the  nominal  case  identified  in  Table  3.4-1. 

The  distance  measurement  accuracies  for  both  atmos¬ 
pheric  monitoring  situations  are  presented  in  Fig.  3.4-3.  For 
single-site  monitoring,  the  Markov  correlation  distance  of 
30  km  (see  Table  3.1-1)  is  traveled  in  240  seconds.  From  the 
upper  curve  in  Fig.  3.4-3,  this  "correlation  time"  corresponds 
to  the  time  at  which  the  position  measurement  accuracy  begins 
to  degrade.  This  degradation  is  necessarily  due  to  a  lack  of 
additional  information  regarding  the  varying  atmospheric  con¬ 
ditions  between  receivers. 

The  sensitivity  of  the  interferometry  analysis  results 
to  unmodeled  atmospheric  property  variations  can  be  reduced  by 
using  more  than  one  refractometer.  Specifically,  multiple 
fixed-site  atmospheric  sensors  can  be  positioned  throughout 
the  survey  area.  Data  collected  from  these  sensors  can  be 
used  to  extrapolate  the  changing  atmospheric  conditions  between 
the  receivers.  An  extension  of  this  multisensor  situation 
involves  dual-site  monitoring  (i.e.,  refractometers  at  both 
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SINGLE-SITE  MONITORING 


TIME  SINCE  SYNCHRONIZATION  (sec) 


Figure  3.4-3  Sensitivity  of  Position  Measurement  Accuracy 

to  Atmospheric  Monitoring 

the  fixed-  and  moving-receiver  locations).  The  results  of 
simulating  this  scenario  are  reflected  in  the  lower  curve  of 
Fig.  3.4-3.  With  atmospheric  monitors  at  both  receivers,  posi¬ 
tion  accuracy  based  on  radio  interferometry  to  GPS  approaches 
75  ppm  rms . 

The  same  effects  are  observed  in  the  acceleration 
measurement  error  results  which  are  illustrated  in  Fig.  3.4-4. 

For  single-sire  monitoring,  the  rms  measurement  error  is  11.0  mgal 
When  the  vehicle  nears  the  correlation  distance  of  the  Markov 
process  (at  240  sec),  acceleration  measurement  accuracy  de¬ 
grades  rapidly.  This  is  expected  due  to  the  lack  of  additional 
information  regarding  varying  atmospheric  conditions  between 
receivers.  Dual-site  monitoring  reduces  the  acceleration  rms 
error  by  nearly  a  factor  of  three. 
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SINGLE-SITE  MONITORING 


TIME  TO  TRAVEL 
CORRELATION  DISTANCE 
AT  460  km/hr 


Figure  3.4-4 


TIME  SINCE  SYNCHRONIZATION  (sec) 

Sensitivity  of  Acceleration  Measurement 
Accuracy  to  Atmospheric  Monitoring 


These  two  simulations  illustrate  "worst-case"  scenarios 
in  which  a  receiver  is  traveling  directly  away  from  an  atmos¬ 
pheric  monitor  positioned  at  a  fixed-site.  As  part  of  the 
preparation  for  analysis  of  radio  interferometry  data  to  be 
recorded  during  a  survey  of  a  particular  area  during  a  particu¬ 
lar  time  of  year,  a  statistical  model  can  be  formulated  to 
describe  accurately  the  varying  atmospheric  conditions  between 
receivers.  For  example,  the  parameters  of  a  Markov  process 
(correlation  distance  and  rms  noise)  can  be  fitted  using  pre¬ 
viously  collected  data  to  portray  the  conditions  expected  at 
the  particular  test  site  during  the  relevant  time-o£-year . 

In  practice,  the  vehicle  trajectory  can  be  carefully  planned 


so  that  the  moving  receiver  is  within  the  correlation  distance 
of  the  fitted  Markov  process  for  sufficiently  long  time  inter¬ 
vals.  For  example,  the  monitor  can  be  positioned  in  the  center 
of  the  area  to  be  surveyed  or  more  than  one  fixed-site  monitor 
can  be  employed  in  the  area.  In  this  situation,  navigation 
accuracies  sufficient  for  the  GGSS  application  would  be  main¬ 
tained  for  longer  times  during  the  course  of  a  survey  mission. 

Preliminary  analyses  indicate  that  the  correlation 
distance  parameter  is  quite  region-sensitive.  Thus,  during 
actual  surveys,  a  single  fixed-site  refractometer  may  not 
be  adequate  to  correct  for  tropospheric  signal  delays  through¬ 
out  the  entire  survey  area.  Furthermore,  as  the  results  in 
Figs.  3.4-3  and  34-4  demonstrate,  single-site  monitoring 
would  not  produce  the  desired  navigation  accuracies.  In  addi¬ 
tion,  dual-site  monitoring  using  microwave  ref ractometers  only 
leads  to  a  factor  of  two  improvement  in  measurement  accuracies 
compared  with  single-site  monitoring.  However,  should  portable 
line-of-sight ,  dual- frequency  radiometers  be  available  for 
survey  use  (see  Section  3.1.1),  vehicle  navigation  accuracies 
can  be  expected  to  improve  by  an  order  of  magnitude  for  dual¬ 
site  monitoring  situations.  Sensitivity  analyses  show  that 
this  improvement  results  in  single-digit  ppm  distance  and  sub- 
milligal  acceleration  measurement  rms  accuracies  for  the  air¬ 
borne  application. 

Sensitivity  of  Measurement  Accuracy  to  Timing  Drift 
Compensation  -  To  compare  the  effects  of  clock  jitter  and  loss 
of  synchronization  between  two  atomic  clocks  on  measurement 
accuracy,  two  different  models  were  employed  in  the  analysis 
(see  Table  3.3-1).  A  random  ramp  process  was  included  to 
model  the  loss  of  synchronization.  The  effect  of  clock  jitter 
on  vehicle  navigation  accuracy  was  studied  by  setting  the  ramp 
slope  parameter  to  zero.  This  latter  case  simulates  the  antici 
pated  post-survey  data- reduction  situation  which  would  occur 
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following  an  accurate  timing  drift  correction.  Additional 
modeling  details  are  discussed  in  Section  3.3  and  summarized 
in  Table  3.3-1;  a  description  of  the  nominal  analysis  scenario 
is  presented  in  Table  3.4-1. 

Although  typical  survey  missions  require  six  to  eight 
hours,  the  effect  of  timing  errors  on  navigation  accuracy  is 
observable  after  only  five  minutes.  The  rms  error  associated 
with  measuring  vehicle  acceleration  is  illustrated  in  Fig.  3.4-5 
for  the  two  situations  described  above.  The  ability  to  measure 
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Figure  3.4-5  Sensitivity  of  Acceleration  Measurement 

Accuracy  to  Adequacy  of  Timing 
Drift  Compensation 
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accurately  acceleration  is  degraded  by  clock  drift  even  when 
atomic  frequency  standards  are  used  and  synchronized  at  the 
start  of  a  survey  mission.  After  five  minutes,  the  accuracy 
of  the  uncorrected  measurements  is  observed  to  be  30  percent 
worse  than  the  case  in  which  the  relative  timing  drift  is 
adequately  compensated.  Thus,  to  ensure  accurate  velocity 
and  acceleration  measurements  using  radio  interferometry,  a 
correction  must  be  applied  to  the  data  to  compensate  for  loss 
of  synchronization. 


Figure  3.4-5  also  illustrates  the  minimal  effect  of 
clock  noise  on  acceleration  uncertainty  after  the  recorded  data 
are  corrected  for  relative  timing  drift.  The  error  due  to 
clock  jitter  ultimately  determines  the  maximum  time  interval 
between  synchronizations.  As  noted  earlier,  survey  operations 
can  last  as  long  as  eight  hours.  During  this  time,  clock  noise 
can  be  expected  to  increase  to  approximately  three  nanoseconds. 
Such  a  timing  error  contributes  a  one  percent  error  in  acceler¬ 
ation  measurement  accuracy  (0.04  mgal  rms  error  at  450  km/hr). 
Accordingly,  an  eight-hour  resynchronization  schedule  is  appro¬ 
priate  for  ensuring  the  negligible  contribution  of  clock  error 
to  navigation  error. 

Sensitivity  of  Measurement  Accuracy  to  Signal  Eleva¬ 
tion  Angle  -  The  sensitivity  of  measurement  accuracy  to  the 
elevation  angle  of  the  satellite-transmitted  signal  was  also 
analyzed.  Signals  from  satellites  near  the  horizon  are  subject 
to  high  refractivity  errors  while  interferometry- related  errors 
increase  whenever  a  tracked  satellite  approaches  the  zenith 
and  the  vehicle  trajectory  becomes  more  parallel  to  the  GPS 
signal  wavefront.  The  effect  of  these  two  extremes  was  studied 
to  determine  the  satellite  "drop"  angle,  i.e.,  the  angle  at 
which  greater  vertical  channel  measurement  accuracy  can  be 
obtained  using  altimeter  rather  than  GPS  information.  Note 
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that  this  drop  angle  is  not  a  priori  GDOP-dependent ,  although 
many  of  the  same  geometric  arguments  apply  (see  Section  3.3.2) 
The  nominal  airborne  scenario  described  in  Table  3.4-1  was 
used  for  this  study. 


The  analysis  results  are  illustrated  in  Fig.  3.4-6. 


As  the  signal  elevation  angle,  0,  increases  beyond  75°,  the 
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acceleration  error  increases  as  1/cos  0 ,  as  predicted  by  the 
discussion  in  Section  3.3.  However,  between  30  and  75  deg  the 
acceleration  rms  error  changes  by  only  0.02  mgal  and  thus,  the 
effect  of  signal  elevation  angle  is  inconsequential  for  the 
practical  case  which  considers  effective  selection  of  avail¬ 
able  satellites. 


3.4.3  Conclusions 


TASC's  analysis  of  moving-receiver  radio  interferom- 


rms  and  70  ppm  rms  are  currently  feasible  for  surface  and  air¬ 
borne  surveys,  respectively.  These  accuracy  levels  exceed 
those  believed  to  be  necessary  for  GGSS  testing.  In  addition, 
acceleration  can  be  measured  to  sub-milligal  accuracy  for 
vehicle  speeds  less  than  100  km/hr,  as  shown  in  Fig.  3.4-2. 

High-accuracy  radio  interferometry-derived  measurements 
will  depend  on  the  manner  in  which  tropospheric  refractivity 
is  measured  and  the  accuracy  with  which  observed  propagation 
delay  is  treated  in  the  system.  Several  refractivity  monitors, 
including  the  class  of  experimental  radiometers,  were  assessed 
and  incorporated  into  simulations  which  replicated  single-  and 
dual-site  atmospheric  monitoring  scenarios.  Not  surprisingly, 
refractivity  monitoring  at  both  receiver  locations  yields  the 
more  accurate  measurements  (compared  with  single-site  monitor¬ 
ing);  dual-site  monitoring  offers  an  improvement  in  navigation 
accuracy  by  a  factor  of  two  or  three.  Thus,  for  certain  survey 
applications  (e.g.,  surveys  local  to  the  monitoring  site), 
detailed  modeling  of  tropospheric  refractivity  errors  may  be 
preferable  to  dual-site  recordings. 

Accurate  inter ferometrically  derived  acceleration 
measurements  require  precise  and  synchronized  clocks  to  main¬ 
tain  data  integrity.  The  analysis  incorporated  error  models 
appropriate  for  atomic  clock  standards  and  examined  the  effects 
of  clock  jitter  and  relative  drift  on  navigation  accuracies. 

The  study  results  indicate  a  need  to  synchronize  the  receiver 
clocks  prior  to  survey  operations  and  every  eight  hours  there¬ 
after,  as  appropriate.  In  addition,  the  recorded  data  must 
be  compensated  for  the  loss  of  clock  synchronization  prior 
to  post-survey  data  reduction. 

The  stud}  also  investigated  available  GPS  receivers 
to  determine  whic  i  is  best  suited  for  moving- receiver  radio 


interferometry  use.  Receivers  were  compared  on  the  basis  of 
hardware  design  and  signal  quantities  measured.  TASC ' s  evalu¬ 
ation  indicates  that  the  Texas  Instruments  TI  4100  GPS  receiver 
is  uniquely  qualified  to  support  applications  of  moving- receiver 
radio  interferometry.  This  result  confirms  the  wisdom  of  the 
current  baseline  GGSS  configuration  which  uses  this  particular 
receiver . 


Based  on  these  and  other  published  results  (Ref.  28), 
the  inter-receiver  distance  and  acceleration  measurement  accur¬ 
acy  levels  presented  herein  can  be  maintained  under  actual 
survey  conditions.  Thus,  it  is  reasonable  to  conclude  that 
moving-receiver  radio  interferometry  offers  an  attractive  backup 
to  the  nominal  airborne  GGSS  positioning  system  (i.e.,  GPS 
pseudoranging)  and  furthermore,  is  a  low-cost,  low-risk  alterna¬ 
tive  to  the  fifth-wheel  mechanization  planned  for  the  GGSS 
surface  vehicle.  The  results  presented  herein  demonstrate  the 
ability  of  radio  interferometry  to  provide  the  stringent  posi¬ 
tioning  accuracies  which  may  be  required  for  surface  testing. 

In  addition,  for  surface  vehicle  surveys,  the  acceleration 
measurements  can  be  used  to  provide  an  independent  measure  of 
gravity  (see  below). 


3.5  TECHNICAL  RISK  AND  RECOMMENDATIONS 
3.5.1  Technical  Risk 

Many  of  the  uncertainties  associated  with  moving- 
receiver  radio  interferometry  have  been  identified  and  quanti¬ 
fied.  These  analysis  results  are  presented  in  the  foregoing 


sections. 


Hardware  considerations  indicate  that  integration  of 
moving-receiver  interferometry  with  the  GGSS  can  be  accomplished 
with  a  minimal  modification  effort  since  the  currently  planned 
GPS  receiver  ard  GGSS  data  acquisition/monitoring  system  are 
adequate  for  p  ovic ing  the  required  phase  measurements.  The 
only  other  additions  required  consist  of  a  second  GPS  receiver, 
atomic  clock,  and  data  recorder  at  a  convenient  fixed  site, 
and  at  least  one  refractivity  monitor  at  either  the  fixed  site 
or  aboard  the  moving  receiver.  Thus,  radio  interferometry  em¬ 
ploying  a  mobile  GPS  receiver  promises  to  be  an  especially 
significant  new  approach  for  providing  very  accurate  navigation 
and  geodetic  information,  either  in  conjuction  with  the  GGSS 
or  as  a  stand-alone  sensor. 

3.5.2  Recommendations 

TASC  has  studied  the  feasiblity  of  using  moving- receiver 
radio  interferometry  as  a  positioning  backup  for  GGSS  surface 
testing  and  survey  operations  and  determined  that  the  approach 
holds  considerable  promise.  These  results  should  be  furthered 
by  augmenting  the  GGSS  surface  test  plan  to  accommodate  this 
capability.  The  1  )gical  next  step  is  a  demonstration  of  inter¬ 
ferometric  acceleration  and  gravity  measurements  in  surface 
vehicle  surveys. 

The  accuracy  potential  for  recovering  the  gravity 
vector  via  multiple  sensor  measurements  bears  further  scrutiny. 

In  particular,  incorporating  inter ferometrical ly  derived  posi¬ 
tion  and  velocity  measurements  with  other  GGSS  aides  (such  as 
a  vertical  accelerometer  or  gravimeter)  may  reduce  airborne 
acceleration  uncertainties  to  sub-milligal  levels  (see  Chap¬ 
ter  4).  Identification  of  the  gravity  field  wavelengths  most 
observable  via  radio  interferometry  is  an  appropriate  extension 
of  this  investigation. 


Finally,  the  status  of  atmospheric  refractivity  moni¬ 
toring  technology  should  continue  to  be  followed  closely.  The 
evaluation  of  rugged,  portable  dual - frequency  microwave  radiom¬ 
eters  will  further  increase  the  capability  of  radio  interferometry 
to  provide  very  accurate  measurements  of  navigation  quantities. 


The  incorporation  of  a  gravimeter  into  the  airborne 
GGSS  as  currently  configured  will  improve  the  quality  of  sur¬ 
face  gravity  disturbance  vector  recovery.  A  gravimeter  is 
complementary  to  the  GGSS  since  it  measures  the  longer  wave¬ 
lengths  of  the  gravity  field,  while  the  GGSS  observes  the 
shorter  wavelengths.  Current  positioning  techniques,  such  as 
those  associated  with  conventional  GPS  pseudoranging,  are 
capable  of  supporting  a  significant  reduction  in  survey  errors 
over  the  entire  gravity  field  spectrum.  A  substantial  improve¬ 
ment  in  navigation  accuracy  (horizontal  velocity  and  vertical 
position,  in  particular),  such  as  may  be  attainable  using  radio 
interferometric  techniques  (see  Chapter  3),  yields  a  reduction 
in  gravity  recovery  errors  over  the  shorter  wavelengths  of  the 
spectrum  (less  than  500  km).  The  associated  velocity  accuracy 
is  necessary  to  overcome  the  dominant  Eotvos  correction  error. 


4.1  SYSTEM  DESCRIPTION  AND  MECHANIZATION  DETAILS 

As  noted  in  Section  2.1,  the  GGSS  mechanization  in¬ 
cludes  three  gravity  gradiometers  mounted  on  a  three-axis 
inertially  stabilized  platform.  Control  updates  based  on 
CPS-determined  position  are  periodically  provided  to  the  plat¬ 
form  to  obtain  adequate  orientation  and  positioning.  The  GGSS 
platform  was  designed  (for  equipment  growth  potential)  to 
accommodate  the  mounting  of  two  Bell  BGM-4  gravimeters  to  the 
pitch  girabal .  Note  that  the  Gravity  Sensor  System,  which  was 
developed  for  the  U.S.  Navy  and  is  the  basis  for  the  GGSS  de¬ 
sign,  includes  two  Bell  BGM-3  gravimeters  In  its  implementation 
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aboard  Che  USNS  Vanguard.  i.n  fact,  for  the  GGSS ,  the  only  changes 
in  the  basic  mechanization  to  support  the  addition  of  the  gravi¬ 
meters  involve  interfacing  between  the  gravimeters  and  the  plat¬ 
form,  data  monitoring/recording  system,  and  onboard  computer. 

A  LaCoste  and  Romberg  air-sea  gravimeter  has  been 
used  as  the  gravity  sensor  in  several  test  and  operational 
airborne  surveys.  In  June  1981,  the  Naval  Research  Laboratory 
tested  an  airborne  gravimetry  system  aboard  a  P3-A  Orion  fixed- 
wing  aircraft  (Ref.  29).  Since  the  start  of  this  decade, 

Carson  Geoscience  Company  (formerly  Carson  Helicopters,  Inc.) 
has  been  conducting  airborne  gravity  surveys  aboard  a  Sikorski 
S61  helicopter  to  aid  oil  and  gas  exploration  (Refs.  30  and  31). 
The  choice  of  this  particular  gravimeter  is  based  primarily  on 
availability  and  cost.  A  brief  description  of  this  gravimeter 
and  two  other  types  which  currently  offer  potential  for  accur¬ 
ately  measuring  gravity  on  a  moving  base  follows. 


The  LaCoste  and  Romberg  gravimeter  uses  a  spring  to 


support  a  pivoted  beam  containing  a  mass  at  its  free  end.  The 
elements  of  the  meter  are  designed  so  that  gravitational  force 
on  the  mass  is  balanced  by  a  "zero- length"  spring.  Such  a 
spring  implies  that  the  tension  is  proportional  to  the  length 
of  the  spring,  and  theoretically,  if  all  external  forces  were 
removed  the  spring  would  collapse  to  zero  length.  In  practice, 
this  gravity  meter  is  used  as  a  null  instrument.  An  adjustable 
second  spring  is  used  to  restore  the  pivoting  beam  to  its 
original  horizontal  position.  The  external  acceleration  can 
then  be  determined  by  measuring  this  resultant  restoring  force 
and  scaling  appropriately.  The  moving  system  of  the  Worden 
gravimeter  is  similar  to  the  LaCoste  and  Romberg  design;  how¬ 
ever,  the  Worden  instrument  is  small,  lightweight,  and  contains 
an  automatic  temperature-compensation  scheme.  Further  details 
of  these  two  gravimeters  are  presented  in  Ref.  32. 


gyrostabilized  platform.  The  sensor  requires  no  cross-coupling 

\ 

correction  and  represents  a  significant  improvement  over  the 
conventional  spring- type  gravimeters.  A  comparison  of  gravity 
data  collected  with  both  the  BGM-3  and  a  conventional- type 
gravimeter  yielded  an  rms  discrepancy  between  free-air  gravity 
anomaly  values  at  ship-track  crossings  of  0.38  mgal  for  the 
BGM-3  and  1.35  mgal  for  the  conventional  type  (Ref.  33).  Based 
on  these  and  other  encouraging  test  results,  the  BGM-3  gravi¬ 
meter  is  now  routinely  used  by  NAVOCEANO  for  gravity  surveys 
over  the  ocean  for  DMA,  and  is  also  used  to  provide  redundant 


gravity  data  in  support  of  ongoing  tests  of  new  equipment.  In 
addition,  a  special  version  of  the  the  BGM-3  is  included  in 
the  Gravimeter  Module  Assembly  (GMA)  of  the  U.S.  Navy's  Gravity 
Sensor  System  (GSS).  Two  GMAs  will  be  mounted  on  the  gravity 
sensor  platform  of  each  GSS.  A  prototype  of  this  version  of 
the  GSS  (Advanced  Development  Model  II)  is  currently  undergoing 
laboratory  testing  at  Bell  Aerospace.  However,  additional 
research  is  needed  prior  to  incorporating  this  more  accurate 
(and  much  more  expensive)  instrument  in  roui ine  airborne  survey 


applications . 


A. 2  ANALYSIS  APPROACH  AND  RESULTS 


A. 2. 1 


Airborne  gravimeter  surveys  are  analyzed  using  TASC ' s 
Multisensor  Survey  Simulation  (MSS)  capability  (Ref.  34). 
This  analysis  represents  the  gravity  field  as  a  two-dimensional 
random  field  process  with  known  correlation  characteristics , 
models  sensor  measurement  errors  as  stationary  gaussian  proc¬ 
esses  independent  of  the  gravity  field,  and  accounts  for  all 


4-3 


significant  geometrical  and  mechanization  effects.  Multi¬ 
dimensional  Wiener  smoothing  is  used  to  derive  the  average 
power  spectral  density  of  post-survey  gravity  residuals.  Since 
the  Wiener  smoother  is  optimal  in  the  sense  of  minimizing  the 
mean-square  survey  error  residuals,  the  MSS  outputs  represent 
the  best  possible  use  of  the  survey  data. 

An  overview  of  the  multisensor  survey  analysis  method¬ 
ology  is  presented  in  Fig.  4.2-1.  The  output  is  obtained  by 
applying  the  analysis  techniques  described  in  Ref.  34.  The 
results  are  the  statistics  of  the  post-survey  residual  gravity 
errors.  The  following  quantites  are  required  inputs  to  the 
analysis  process: 

•  Statistical  model  for  the  unsurveyed 
anomalous  field 

•  Gravimeter  and  GGSS  error  models  appro¬ 
priate  for  an  airborne  application 

•  Characteristics  of  the  survey. 


For  this  study,  a  statistical  model  for  gravity  in 
the  North  Texas  GGSS  test  area  was  employed.  The  model  is  an 
Attenuated  White  Noise  (AWN)  gravity  field  model  which  was 
specially  tuned  to  the  physiography  of  the  GGSS  test  area. 
Additional  details  of  this  model  are  presented  in  Ref.  35. 

4.2.2  Survey  Geometry  and  Measurement  ilrror  Models 

Survey  geometry  refers  to  the  position  of  the  measure¬ 
ments  and  the  orientation  of  the  sensors  with  respect  to  the 
earth.  Measurement  error  models,  discussed  in  more  detail 
below,  are  used  to  characterize  the  applicable  error  sources. 
The  airborne  gravimeter  and  GGSS  error  models  were  derived  by 
appropriately  extending  the  shipborne  gravimetry  and  airborne 
gradiometry  models  documented  in  Ref.  34.  The  discussion  below 
follows  the  development  and  conventions  used  in  Refs.  34  and  36. 

Airborne  GGSS  Survey  -  Figure  4.2-2  identifies  the 
nominal  GGSS  airborne  survey  parameters  which  were  used  to 
estimate  the  statistics  of  the  residual  gravity  errors.  The 
cross-track  spacing  is  x ^  anc*  the  along-track  spacing  is 

=  Vtg  (4.2-1) 

where  V  is  the  nominal  survey  speed  and  t  (10  sec)  is  the 
time  interval  between  successive  samples. 

As  shown  in  Fig.  4.2-3,  the  three  gravity  gradiometer 
instruments  are  symmetrically  distributed  about  the  vertical 
axis,  with  each  gradiometer  inclined  at  the  same  "umbrella" 
angle  a.  Figure  4.2-3  also  identifies  the  u,v,w,  instrument 
axes  and  local-level  x,y,z  platform  axes.  Each  gradiometer 
measures  two  elements  of  the  gravity  gradient  tensor  in  the 
instrument  u,v,w  coordinate  frame,  i.e.,  the  cross  gradient 
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Tjv  and  the  difference  of  the  two  in-line  gradients,  Tuu  -  Tvv> 
Furthermore,  each  of  these  two  gradiometer  outputs  is  a  certain 
linear  combination  of  the  gradients  of  the  potential  given  by 

"w  r -i/a  1/4  o 

T1?  i  -1/3  0  1/3 

i 

T>,  1/4  -1/4  0 

: , ,  -1/3  0  1/3 

r ,  '  o  o  o 

Tj2  1/6  -1/2  1/3 

(4.2-2) 

where  represents  the  j  output  of  the  i  gradiometer 

<j  =  !'2'  *"d  1  =  3,2,3)  and  T^.  Tyy,  Tzz ,  Txy,  Txz ,  and  Tyz 
are  elements  of  the  gravity  gradient  tensor  in  an  east  (x), 
north  (y),  up  (z)  reference  frame. 

The  vector  transfer  function,  F,  from  anomalous  surface 
potential  to  the  measured  quantities  can  be  represented  by 

F(s)  =  BH(s)  (4.2-3) 

where  s  is  the  vector  of  spatial  frequencies  (s^  and  S2  are 
measured  in  the  east  and  north  directions,  respectively);  B 
is  the  6x6  matrix  of  Eq.  4.2-2;  and  H,  the  vector  transfer 
function  from  anomalous  potential  at  height  h  is  given  by 


J5/6  VT/2  yfZ/6 

-yfT/3  -V7/6  W  6  Tyy 

S5/6  -J2/2  ft/6  Tzz 

VT/3  -ft/6  -ft/6  Txy 

-V7/3  0  -ft/3  Txz 

0  ft/3  0  T 

J  L  yz 


H(s)  =  4n2e'27ish 


where  s 


s22)1/2 


-sls2 

-  is^s 
~^S2S 


(4.2-4) 


The  spectral  density  of  the  measurement  errors  is 
derived  next.  Analysis  by  Bell  Aerospace  of  test  data  from 
the  gradometers  which  make  up  the  GGSS  yielded  the  following 
analytic  form  for  the  self-noise  spectra: 

SN,N(f)  =  r/f1'6  +  w  (4.2-5) 

2 

where  f  is  measured  in  Hz,  r  is  the  red-noise  level  in  E  x  Hz, 

2 

and  w  is  the  white-noise  floor  in  E  /Hz. 


Temporal  frequency,  f,  is  transfered  into  spatial 
frequency,  s-^,  through  the  mapping 


f  = 


(4.2-6) 


On  different  data  tracks,  measurement  errors  are  taken  to  be 
independent  to  give  a  diagonal  error  spectral  density  matrix, 
$£  jr(s)  ,  with  entries 


*  p  p  ( I*) 
K  ’  K 


T  2V 


rk/V‘ 


+  w, 


(4.2-7) 


for  -l/(2t^)  <  <  l/(2x^)  and  k  =  1,2, ••■,6.  The  values  of 

and  w^  which  were  used  in  the  simulation  are  presented  in  Table  4.2 
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TABLE  4.2-1 

RED-NOISE  AND  WHITE-NOISE  PARAMETERS 


QUANTITY 

DOUBLE -SIDED 
PSD 

SINGLE-SIDED 

PSD 

Red-Noise 

Level 

1.7  x  10"4  E2Hz 

1.0  x  10'3  E2  rad/sec 

White-Noise 

Floor 

_ 

55  E2/Hz 

17.5  E2/( rad/sec) 

Note:  The  single-sided  PSD,  F  ,  is  defined  for  frequencies 

greater  than  zero.  The  double-sided  PSD,  ,  is 

defined  for  both  positive  and  negative  frequencies, 
but  is  an  even  function.  For  positive  frequencies, 
the  relation  between  Fg  and  F^  is 


/  0  ( u) )  du> 

0 

or  equivalently,  Fg  =  Fp/n,  where  f  and  tu  are 
the  frequency  in  Hz  and  rad/sec,  respectively. 

Airborne  Gravimeter  Survey  -  The  geometry  of  an  air¬ 
borne  gravimeter  survey  was  selected  to  be  identical  to  that 
of  the  bidirectional  GGSS  survey  illustrated  in  Fig.  4.2-2. 
The  cross-track  spacing,  t2»  is  equal  to  5  km  and  the  along- 
track  spacing  is 

=  V  At  (4 . 2-8) 

where  At,  the  time  interval  between  successive  samples,  was 
chosen  to  be  12  sec.  This  sampling  time  corresponds  to  a  1  -  km 
data  spacing  in  the  along-track  direction. 


=  (1 /n)  f  0 ( f )df 
0 
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The  transfer  function,  F(s),  from  surface  anomalous 
potential  to  the  measured  gravity  quantites  was  derived  in 
Ref.  34  for  the  shipborne  gravimeter  survey  scenario  and  is 
identical  in  form  to  the  aircraft  case  being  considered  here. 
The  expression  for  1 (^)  is 

F(s)  =  (2ns  -  2/R)  +  2Qcos0 [ i2ns^  ]  G(s^)  (4.2-9) 

where  the  two  terms  are  the  transfer  functions  from  anomalous 
surface  potential  to  gravity  anomaly  and  east  component  of  the 
gravity  disturbance  vector  at  the  earth's  surface,  respectively, 
and 

R  =  the  earth  radius  (6.378  *  10^  m) 

0  =  the  earth  rotation  rate  (0.7292  x  10  ^  rad/sec) 

0  =  the  reference  latitude 

G(s^)  =  the  navigation  transfer  function  from  sensed 

acceleration  to  indicated  velocity  as  a  function 
of  distance  in  the  east  direction,  given  by 


s1/(2nV) 

G(  s  i  )  =  - - - 7 - 5-  (4.2-10) 

2p(fo/V)Sl  +  i [ s :  - ( f Q/V ) ^ ] 

p  =  the  navigator  damping  coefficient 
f  =  the  Schuler  frequency  (1.971  x  10  ^  Hz). 

The  airborne  survey  measurement  error  model  consists 
of  three  error  source  categories:  1)  instrument  errors,  2)  ver¬ 
tical  position  uncertainty,  and  3)  Eotvos  correction  errors. 
These  three  categories  are  each  discussed  below. 

The  gravimeter  recordings  are  corrupted  by  various 
error  sources,  including  quantization  and  instrument  noise. 


These  effects  are  combined  into  a  single  term,  I,  which  is 
modeled  as  being  independent  from  measurement  point  to  point 
and  track  to  track.  Correlated  error  sources  (e.g.,  errors  in 
the  Eotvos  correction)  are  treated  separately  as  indicated 
below.  An  appropriate  rms  value  for  the  LaCoste  and  Romberg 
gravimeter  in  an  airborne  application  is  =  1.0  mgal. 

The  PSD  of  I  is 

<t> i  -j^s)  =  1  20I^  (A. 2-11) 

Vertical  position  uncertainty  affects  the  gravity 
survey  results  by  corrupting  the  calculated  vertical  gradient 
of  the  field  and  the  aircraft  vertical  acceleration  correction 
as  explained  in  Ref.  29.  The  combined  effect  of  this  error  is 
represented  by  a  single  term,  H,  which  is  modeled  as  being 
independent  from  measuremer *  point  to  point  and  track  to  track. 
An  appropriate  rms  value  for  this  source  in  the  currently  con¬ 
figured  GGSS ,  where  vertical  position  will  be  provided  by  a 
GPS-inertial  or  barometric  altimeter  mechanization  is  =  20  m. 
The  PSD  of  H  is 

,H^— ^  =  xlT2^2go/,R^  °H2  (A. 2-12) 

where  2gQ/R  =  0.308  mgal/m. 

The  Eotvos  correction  error  is  caused  by  the  error  in 
the  estimate  of  the  east  component  of  velocity  obtained  from 
the  inertial  navigator  aboard  the  survey  aircraft.  This  east 
velocity  error  is  attributable  to  the  east  component  of  the 
gravity  disturbance  vector  sensed  by  the  navigator  along  with 
random  gyro  and  accelerometer  errors.  Thus,  if  airborne  gravi¬ 
metric  data  are  viewed  as  gravity  anomaly  measurements,  the 
measurement  errors  are  correlated  with  the  gravity  field. 


I 

* 

► 


However,  ar  equivalent  formulation  can  be  obtained  when 
the  data  are  viewed  as  noisy  measurements  of  linear  combinations 
of  the  gravity  anomaly  and  the  east  component  of  the  gravity 
disturbance  vector  modulated  by  the  navigator  response.  This 
convention  is  adopted  in  Ref.  34  and  is  also  used  here  since, 
in  this  formulation,  measurement  errors  turn  out  to  be  indepen¬ 
dent  of  the  gravity  field  (see  Ref.  34). 

Using  this  formulation,  the  Eotvos  correction  error 
source,  C,  is  modeled  as  a  white-noise  sensed  acceleration 
input  with  constant  spectral  density  q.  The  magnitude  of  q 
is  given  by 

q  =  8nV  fQ  po^  (4.2-13) 

where  a v  is  the  rms  velocity  error  which  is  about  0.2  m/sec 
for  the  currently  configured  GGSS.  Following  the  development 
in  Ref.  34,  the  spectral  density  of  the  uncorrelated  Eotvos 
correction  errors  in  the  interval  -1/(2t^)  <  s-^  <  l/(2x^)  is 

<t>c  c(s)  =  ^fi^cos^9x  2^  I  ^(  si  ^  (4.2-14) 


In  summary,  the  spectral  density  of  the  airborne 
gravimeter  survey  measurement  errors,  p£  £.  is  given  by  the 
sum  of  the  spectral  densities  of  the  instrument-related  errors 
(Eq.  4.2-11)  and  the  errors  due  to  vertical  position  uncertainty 
(Eq.  4.2-12)  and  velocity  uncertainty  (Eq.  4.2-14).  Thus, 


<J,E,E(-)  = 


I(s) 


♦h,h(s) 


c.c 


(  S  ) 


(4.2-15) 


for  -1/(21,)  <  s,  <  l/(2t,)  and  -1/(2t9)  <  s0  <  l/(2t0). 


4.2.3  Results 

This  section  describes  the  improvement  in  survey 
error  which  can  be  anticipated  by  incorporating  a  gravimeter 
into  the  nominal  GGSS  mechanization.  The  results  which  follow 
were  obtained  using  the  multisensor  survey  analysis  and  the 
engineering  error  models  discussed  in  the  previous  section. 
Likewise,  the  nominal  airborne  scenario  and  the  North  Texas 
gravity  field  model  were  also  identified  in  Section  4.2.2. 

Table  4.2-2  contains  the  rms  survey  error  for  three 
particular  survey  configurations:  1)  nominal  GGSS,  2)  nominal 
GGSS  with  gravimeter  aiding,  and  3)  improved  GGSS  with  gravi¬ 
meter  aiding.  The  nominal  GGSS  case  is  the  baseline  configura¬ 
tion  described  in  Chapter  1  (i.e.,  GGSS  with  GPS  aiding).  The 


TABLE  4.2-2 

GRAVIMETER-AIDING  ANALYSIS  RESULTS 


SURVEY 

CONFIGURATION 

POST-SURVEY  RMS  DEFLECTION 

OF  THE  VERTICAL  ERROR 
(arc  sec)  (EACH  AXIS) 

POST-SURVEY  RMS  VERTICAL 
GRAVITY  DISTURBANCE  ERROR 
(mgal) 

ALL 

WAVELENGTHS 

(A) 

A  <  500  km 

ALL 

WAVELENGTHS 

(A) 

A  <  500  km 

Performance 

Goal 

0.20 

0.18 

1.0 

0.90 

0.12 

3.5 

0.75 

Nominal  GGSS 
with  Gravi¬ 
meter  Aiding 

0.13 

0.12 

0.84 

0.75 

Improved  GGSS 
with  Gravi¬ 
meter  Aiding 

0.11 

0.11 

0.71 

0.71 

'"Tor  the  case  in  which  no  tie-point  data  are  available. 
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nominal  GGSS  with  gravimeter  aiding  case  represents  modifying 
the  baseline  GGSS  mechanization  (as  discussed  in  Section  4.1) 
to  accommodate  implementation  of  a  gravimeter  on  the  GGSS  plat¬ 
form.  Finally,  the  improved  GGSS  with  gravimeter  aiding  case 
reflects  the  incorporation  of  very  accurate  navigation  informa¬ 
tion  into  the  nominal  GGSS  with  gravimeter  aiding  mechanization. 

The  post  survey  rms  error  represents  the  contribution 
of  all  known  error  sources,  including  sampling  effects,  down¬ 
ward  continuation,  finite  extent,  and  data  processing  algorithm 
errors.  The  error  allocation  procedure  outlined  in  Ref.  36 
was  used  to  obtain  these  results.  Also  indicated  in  this  table 
for  comparison  purposes,  is  DMA's  overall  performance  goal  for 
estimation  of  the  gravity  disturbance  vector  components  at  the 
earth's  surface.  The  entries  in  Table  4.2-2  are  separated  into 
two  columns  for  post-survey  rms  errors  over  all  wavelengths 
and  for  only  wavelengths  less  than  500  km.  This  separation  is 
consistent  with  DMA's  performance  goal  which  was  imposed  to 
accommodate  manageable  testing  of  the  GGSS  within  its  perform¬ 
ance  limits. 

For  gravity  field  wavelengths  shorter  than  500  km, 
adding  a  gravimeter  to  the  nominal  GGSS  configuration  offers 
no  improvement  to  survey  accuracy.  However,  a  substantial  im¬ 
provement  in  survey  accuracy  is  realized  at  longer  wavelengths. 
In  fact,  DMA's  overall  performance  goal  (1.0  ragal  rms  recovery 
accuracy  over  all  wavelengths)  appears  achievable  with  the 
gravimeter-aided  GGSS  configuration.  This  advantage  attributes 
to  the  gravimeter  being  complementary  to  the  GGSS.  The  gravi¬ 
meter  accurately  measures  the  longer  wavelengths  of  the  gravity 
field  while  the  GGSS  accurately  observes  the  shorter  wavelengths 


-'■Currently,  GGSS  surveys  are  envisioned  to  incorporate  surface 
"tie-point"  gravity  data  as  a  source  of  long-wave  length 
i u formation . 


Thus,  it  is  not  surprising  that  no  improvement  is  obtained  at 
the  shorter  wavelengths  by  adding  a  gravimeter  to  the  GGSS 
mechanization .  The  GGSS  already  recovers  this  part  of  the 
gravity  field  spectrum  accurately  and  a  gravimeter  does  not. 

Sensitivity  analyses  confirm  that  the  vertical  posi¬ 
tion  uncertainty  is  a  significant  contributor  to  airborne 
gravimeter  survey  error.  Furthermore,  the  Eotvos  correction 
error  is  also  an  important  and  distinct  error  source.  The 
fourth  survey  configuration  in  Table  4.2-2  is  included  to 
demonstrate  the  improvement  in  survey  accuracy  which  can  be 
expected  when  these  two  effects  are  reduced  substantially. 

In  this  particular  case,  the  vertical  position  rms  error  was 
reduced  to  below  1.0  cm  and  the  velocity  rms  error  to  2  cm/sec. 
These  improvements  in  navigation  accuracy  yield  about  five 
percent  improvement  in  survey  accuracy  for  gravity  field  wave¬ 
lengths  shorter  than  500  km  and  about  15  percent  improvement 
over  all  wavelengths.  The  case  in  Table  4.2-2  with  which 
these  reductions  are  compared  is  the  nominal  GGSS  with  gravi¬ 
meter  aiding. 

To  attain  the  good  navigation  accuracies  noted  above, 
additional  mechanization  changes  in  the  nominal  GGSS  are  re¬ 
quired.  Discussions  of  the  use  of  high- resolution  radar, 
laser,  and  pressure  altimeters,  radar  terrain  clearance,  and 
multiple  electronic  navigation  tracking  techniques  to  reduce 
the  vertical  position  uncertainty,  are  presented  in  Refs.  29 
and  31.  The  use  of  moving-receiver  radio  interferometry  is 
another  possibility  which  is  especially  promising  since  it 
offers  an  improvement  in  navigation  accuracy  in  all  three 
axes.  Additional  details  of  moving- receiver  radio  interfer¬ 
ometry  are  presented  in  Chapter  3. 
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4.3 


TECHNICAL  RISK  AND  RECOMMENDATIONS 


The'  expected  improvement  in  survey  accuracy  over  all 
gravity  field  wavelengths,  which  would  be  realized  by  incorpo¬ 
rating  a  gravimeter  into  the  nominal  GGSS  configuration,  moti¬ 
vates  continuing  research  in  this  area.  The  primary  technical 
risk  area  relates  to  the  inclusion  of  an  appropriate  gravi¬ 
meter  on  the  GGSS  inertially  stabilized  platform.  The  platform 
was  designed  to  accept  the  Bell  BGM-4  model  gravimeter,  but 
consideration  has  also  been  given  to  adapting  the  Bell  BGM-3 
model . 

In  addition,  technical  risk  is  associated  with  obtain¬ 
ing  the  navigation  accuracies  needed  for  an  improvement  in 
survey  accuracy  over  the  shorter  wavelengths  of  the  gravity 
field.  Navigation  accuracies  of  that  level  are  currently 
possible  only  by  exploiting  state-of-the-art  concepts.  Addi¬ 
tional  analysis  and  demonstrations  of  these  concepts  are 
recommended . 

In  particular,  TASC  recommends  that  a  demonstration  * 
of  radio  interferometry  in  a  moving- receiver  application  be 
conducted  in  the  near  future.  This  demonstration  would  deter¬ 
mine  if  the  navigation  accuracies  predicted  herein  can  be 
obtained  in  a  practical  manner  under  "real-world"  operating 
conditions.  Upon  successful  demonstration  of  the  moving- 
receiver  radio  interferometric  concept,  consideration  should 
be  given  to  incorporating  a  gravimeter  into  the  GGSS  mechani¬ 
zation  to  determine  if,  in  fact,  relief  can  be  obtained  from 
the  burden  of  emplacing  surface  gravity  and  deflection  tie 
points  in  support  of  GGSS  surveys. 
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MASTER  INS  ANALYSIS 


The  advantages  of  a  high-accuracy  inertial  navigation 
system  (master  INS)  used  to  augment  the  airborne  GGSS  were 
examined  as  part  of  the  aided-airborne  GGSS  study.  The  exam¬ 
ination  addressed  both  the  aircraft  navigation  and  gradiometer 
platform  stabilization  requirements.  Several  state-of-the-art 
inertial  systems  which  possess  position  error  rate  specifica¬ 
tions  of  better  than  0.3  nm/hr  CEP  were  considered  as  potential 
GGSS  aides.  A  list  of  these  systems  and  summary  of  their  prac¬ 
ticality  for  GGSS  aiding  are  presented  in  Table  5-1. 

TABLE  5-1 

SUMMARY  OF  MASTER  INS  EVALUATION 


SYSTEM 

COMMENTS 

AIRS 

COSTLY;  EXTENSIVE  EFFORT  REQUIRED  FOR  USE  IN  SURVEY 

AIRCRAFT 

LN-37 ,  CAROUSEL, 

N-73  AND  SKN-2440 

EVEN  WITH  VERTICAL  DEFLECTION  COMPENSATION,  UNAIDED 

INERTIAL  OPERATION  WILL  PROVIDE  NAVIGATION  ACCURACIES 

LESS  THAN  200  m  ONLY  FOR  ABOUT  ONE  HOUR 

SPN/GEANS 

SLIGHTLY  BETTER  THAN  ABOVE ;  ABOUT  1 . 5  HOURS  BEFORE 

NAVIGATION  ERRORS  EXCEED  200  m 

NAS-26 

i 

BASED  ON  LIMITED  VAN  TEST  RESULTS,  ADEQUATE  (UP  TO 

16  HOURS)  OPERATIONS  ARE  THEORETICALLY  POSSIBLE;  FLIGHT 

TEST  RESULTS  YIELDED  LESS  ENCOURAGING  RESULTS;  AIRCRAFT 

WOULD  HAVE  TO  BE  CONFIGURED  WITH  STELLAR  VIEWING  APERTURE; 

CLOUD  COVER  WOULD  LIMIT  OPERATIONS 

1 
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The  key  point  associated  with  using  a  master  INS  is 
that  in  survey  areas  where  the  gravity  disturbances  are  typical 
of  the  worldwide  average,  no  purely  inertial  system  can  meet 
the  GGSS'  100  m  CEP  requirement  for  the  duration  of  a  mission 
(at  least  four  hours)  without  some  form  of  compensation  for 
those  gravity  effects.  This  is  true  even  for  a  "perfect" 
astroinertial  system.  (Note  that  gravity-induced  errors  are 
not  observable  in  the  star  sightings.)  Moreover,  even  by  using 
the  gradiometer  to  perform  real-time  (or  post-survey)  vertical 
deflection  compensation,  the  inherent  accuracy  of  currently 
planned  gradiometer  instruments  is  not  adequate  to  support 
long-term  autonomous  navigation.  Under  the  most  optimistic 
assumptions,  only  about  one  hour  of  sufficiently  accurate 
unaided  inertial  navigation  would  be  possible,  compared  with  5 
to  10  minutes  for  the  presently  planned  GPS-aided  GGSS  con¬ 
figuration.  This  modest  extension  in  operating  time  does  not 
justify  the  extensive  hardware/software  modifications  which 
would  be  required  to  incorporate  the  master  INS  capability. 

Thus,  based  on  the  results  of  this  evaluation  and 
on-going  *  navigation  analysis  of  the  baseline  GGSS  configura¬ 
tion,  use  of  a  high-accuracy  INS  as  a  GGSS  aide  is  unnecessary. 
Additional  support  material  and  technical  details  of  the  asso¬ 
ciated  evaluation  are  presented  in  Ref.  1. 


CONCLUSIONS  AND  RECOMMENDATIONS 


This  report  describes  the  currently  planned  GGSS  im¬ 
plementation  along  with  several  variants  which  were  considered 
to  assure  long-term  success  of  the  GGSS.  In  particular,  a 
moving-receiver  radio  interferometer,  gravimeter,  and  master 
INS  are  addressed.  The  sections  which  follow  summarize  TASC's 
key  conclusions  and  identify  the  areas  most  promising  for 
continued  development. 


6.1  CONCLUSIONS 


•  Conventional  GPS  pseudoranging  techniques 
are  more  than  adequate  for  supporting 
test  and  operational  GGSS  missions  (the 
short  outages  associated  with  the  six- 
plane,  18-satellite  constellation  are 
due  to  poor  geometry  and  are  inconsequen¬ 
tial  to  GGSS  deployment  as  currently 
configured) 

•  GPS  backup  modes  based  on  the  use  of  an 
atomic  clock  and  moderately  accurate 
altimeter  will  provide  adequate  extended 
coverage  and  geometry  appropriate  for 
GGSS  navigation  during  testing 

•  Moving-receiver  radio  interferometry 
offers  a  low-cost,  low-risk  backup  posi¬ 
tion  measurement  capability  especially 
well-suited  for  GGSS  surface  testing 

•  The  GGSS  Program  offers  a  unique  oppor¬ 
tunity  for  demonstrating  the  gravity 
measurement  capability  associated  with 
moving-receiver  radio  interferometry 


•  For  short-wavelength  gravity  recovery, 
substantial  improvements  in  navigation 
accuracy  (such  as  would  be  provided  by 
radio  interferometry  to  GPS)  are  neces¬ 
sary  for  gravimeter  aiding  to  yield 
better  performance  than  that  possible 
from  a  nominal  GGSS  survey,  as  well  as  a 
reduction  in  required  surface  tie-point 
data 

•  State-of-the-art  master  INSs  offer  dimin¬ 
ished  incremental  returns  in  autonomous 
navigation  compared  with  the  baseline 
GPS-aided  GGSS  configuration. 


6 . 2  RECOMMENDATIONS 


TASC  recommends  the  following  next  steps: 


•  Specify  appropriate  testing  procedures 

to  assure  confidence  in  the  GPS  Navigator 
during  GGSS  testing  (see  Section  2.1) 

•  Incorporate  moving-receiver  radio  inter¬ 
ferometry  into  the  GGSS  Program  as  a 
backup  position  measurement  capability 
during  surface  testing  (see  Section  3.4.2); 
to  minimize  program  delays;  at  the  present 
time,  configure  only  the  base  station  with 
the  necessary  hardware  (GPS  receiver/recorder, 
atomic  clock,  weather  sensor)  -  the  moving 
vehicle  can  be  augmented  later,  if  necessary 

•  Augment  the  GGSS  test  program  to  accom¬ 
modate  demonstrations  of  gravity  measure¬ 
ment  using  radio  interferometric  principles 
(see  Section  3.4.2) 

•  Support  the  development  of  weather  sen¬ 
sors  which  will  very  accurately  determine 
the  spatial  variations  of  tropospheric 
refractivity  (see  Section  3.4.2) 

•  Perform  application-oriented  analysis  of 
moving-base  gravimetric  surveys  (both 
theoretical  analysis  and  data  processing)  - 
see  Section  4.2.2. 
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APPENDIX  A 

GPS  DESCRIPTION  AND  ANALYSIS  RESULTS  DISCUSSION 


This  appendix  contains  a  detailed  discussion  of  the 
Global  Positioning  System  (GPS).  Section  A.l  describes  the 
orientation  of  the  satellites  which  make  up  the  space  segment 
and  defines  the  attributes  of  the  transmitted  signals  which 
will  be  received  by  the  user  segment.  The  conversion  of  these 
signals  into  positioning  and  navigation  information  is  explained 
and  the  significance  of  satellite  geometry  on  navigation  accur¬ 
acy  is  quantified. 

Section  A. 2  treats  the  analytical  problem  of  modeling 
GPS  satellite  orbits.  Basic  orbital  relationships  and  param¬ 
eters  are  established  in  terms  of  idealized  elliptical  orbits, 
and  perturbation  theory  is  introduced  to  account  for  the  effect 
of  the  ellipsoidal  shape  of  the  earth  on  satellite  orbits. 
Common  features  of  GPS  satellite  orbits  are  examined,  and  the 
problem  of  simulating  a  constellation  of  satellites  is  considered 

Section  A. 3  presents  GPS  analysis  results  involving 
positioning  and  navigation  accuracies  which  can  be  expected 
when  the  complete  18-satellite  constellation  is  operational. 

Note  that  the  GPS  operational  schedule  corresponds  to  the 
currently  planned  timeframe  for  deployment  of  the  GGSS  for 
operational  surveying.  Predicted  satellite  visibility  and 
geometry  are  emphasized  throughout  the  discussion  of  the  an¬ 
alysis  results. 


A.l  SYSTEM  DESCRIPTION 

The  Space  Segment  -  Current  plans  call  for  the  deploy¬ 
ment  of  18  primary  GPS  satellites  in  six  orbital  planes.  Each 


NUMBER  OF  AVAILABLE  SATELLITES 


plane  will  be  inclined  at  55  deg  to  the  equator.  In  addition, 
three  spares  will  be  placed  in  orbit  to  insure  a  high  degree 
of  system  availability.  The  current  schedule  for  deployment 
is  provided  in  Fig.  A. 1-1  and  is  based  on  material  presented 
in  Ref.  8.  The  current  airborne  GGSS  testing  schedule  is  also 
indicated  on  this  figure.  As  is  clear  from  Fig.  A. 1-1,  back-up 
navigation  modes  must  be  considered  to  assure  efficient  and 
successful  GGSS  testing. 

The  current  test  configuration  consists  of  six  "system 
validation"  or  Block  I  satellites.  These  satellites  are  in 
orbits  inclined  at  63  deg  and  are  positioned  to  provide  good 
geometry  for  a  user  at  the  Yuma  Proving  Ground.  Table  A.  1-1 
summarizes  the  Block  I  GPS  orbital  parameters;  Section  A. 2 
defines  the  associated  orbital  parameters  and  explains  their 
significance . 
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Figure  A. 1-1  GPS  Satellite  Deployment  Schedule 
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TABLE  A. 1-1 

BLOCK  I  GPS  ORBITAL  PARAMETERS 


T- 53  71 


POSITION 

NUMBER 

RIGHT  ASCENSION 
OF  NODE 
(deg  ±  2  deg) 

LONGITUDE  OF 
NODE 

(deg  E  t  2  deg) 

TIME  OF  NODE 
(hr : min : sec : GMT 
i  16  min ) 

r 

INCLINATION 
( deg  t  1  ) 

... 

ANOMALY 
(deg  i  2  deg) 

1 

202 

227.0 

3:40:16 

63 

53 

2 

202 

248.5 

2:19:25 

63 

96 

3 

202 

269.5 

0:51:18 

63 

138 

4 

82 

130.5 

2:11:26 

63 

97 

5 

82 

152.0 

0:40:51 

63 

1  AO 

6 

82 

353.0 

11:15:21 

63 

182 

Note:  Values  apply  to  time  at  first  ascending  node  referred  to  1  July  1980. 


Transition  from  the  Block  I  to  the  operational  (Block 
II)  configuration  involves  considering  an  appropriate  build-up 
strategy.  This  strategy  must  take  into  account  ongoing  testing 
requirements  at  Yuma  and  the  launch  windows  associated  with  the 
planned  missions  of  the  space  shuttle  which  currently  is  desig¬ 
nated  as  the  launch  platform  for  these  satellites.  One  such 
strategy  is  described  in  Ref.  37  and  illustrated  in  Fig.  A. 1-2. 
The  rephasing  of  the  Block  I  satellites  into  their  Block  II 
positions  shown  in  Fig.  A. 1-2  is  based  on  this  particular 
strategy . 


Table  A. 1-2  summarizes  the  Block  II  orbital  parameters. 
Three  satellites  will  be  deployed  in  each  of  six  orbital  planes 
which  will  be  equally  spaced  60  deg  apart  in  longitude.  The 
satellites  in  each  plane  will  also  be  uniformly  spaced  120  deg 
apart.  The  relative  phasing  of  satellites  from  one  plane  to 
the  next  will  be  40  deg  (i.e.,  a  satellite  in  one  plane  will 
have  another  satellite  40  deg  "ahead"  or  North  of  it  in  the 
plane  directly  to  the  east,  as  shown  in  Fig.  A. 1-2.  The  bene¬ 
fits  associated  with  this  six-plane  arrangement  of  satellites, 
in  terms  of  continuous  intervals  of  geographic  coverage  and 
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PLANE 


•  NOMINAL  LOCATION  AND  POSITION  NUMBER  OF  BLOCK  II  SATELLITES 
O  NOMINAL  LOCATION  AND  POSITION  NUMBER  OF  BLOCK  I  SATELLITES 
6^  PROPOSED  MOVE  OF  A  BLOCK  I  SATELLITE  TO  A  BLOCK  II  POSITION 


Figure  A. 1-2  A  Strategy  for  Buildup  of  Operational 

GPS  Configuration  (Ref.  37) 

geometry  of  the  satellites  in  view  relative  to  the  users,  are 
discussed  in  Refs.  38  and  39. 

The  satellites  will  be  positioned  in  near-circular 
orbits  at  altitudes  of  approximately  20,183  km.  This  altitude 
was  selected  to  provide  a  period  exactly  one-half  synchronous 
(i.e.,  11  hr  57  min  58.3  sec).  Thus,  each  satellite  will  pass 
over  the  same  point  on  the  earth  every  23  hr  55  min  56.6  sec 
(a  sidereal  day),  or  in  other  words  the  satellites  will  pass 
over  the  same  location  about  four  minutes  earlier  each  day. 

The  User  Segment  -  Each  GPS  satellite  continuously 
broadcasts  a  unique  navigation  message  which  can  be  used  by 
the  user  equipment  to  determine  the  navigation  solution.  The 
message  consists  of  transmitted  signals  characterized  by  two 
components : 


TABLE  A. 1-2 

BLOCK  II  GPS  ORBITAL  PARAMETERS 


ORBITAL 

PLANE 

NUMBER 

SATELLITE 

POSITION 

NUMBER 

RELATIVE 

ANOMALY 

(deg) 

RELATIVE 

RIGHT 

ASCENSION 

(deg) 

RIGHT  ASCENSION 
OF  ASCENDING 
NODE  (deg)* 

— 

LONGITUDE  OF 
ASCENDING 

NODE  (deg) 

1 

0 

0,180 

1 

2 

120 

0 

30 

240,60 

3 

240 

300,120 

4 

40 

260,80 

2 

5 

160 

60 

90 

320,140 

6 

280 

20,200 

7 

80 

340,160 

3 

8 

200 

120 

150 

40,220 

9 

.  .  _  _ 

320 

100,280 

10 

120 

60,240 

4 

11 

240 

180 

210 

120,300 

12 

0 

• 

180,0 

13 

160 

140,320 

5 

14 

280 

240 

270 

200,20 

15 

40 

80,260 

16 

200 

220,40 

6 

17 

320 

300 

330 

280,100 

18 

80 

160,340 

1 

19 

30 

0 

30 

95,15 

5 

20 

310 

240 

270 

215,35 

3 

21 

170 

120 

150 

25,105 

•Referenced  to  astronomical  coordinates  of  1950.0  as  of  0  hr  0  min  GMT  on 
1  July  1985  and  regressing  at  -0.04009  deg/day. 
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(1)  A  coarse/acquisition  (C/A)  signal  and  a 
precise  (P)  signal  centered  about  the  Link  1 
L-band  (L^)  carrier  frequency 

(2)  Either  the  C/A  or  P  signal  centered  about 
the  Link  2  L-band  (L2)  carrier  frequency. 

The  carrier  frequencies  were  selected  to  be  integral  multiples  of 
the  basic  satellite  clock  frequency  of  10.23  MHz.  In  particular, 

L^  =  1575.42  MHz  (19-cm  wavelength)  =  154  x  10.23  MHz 

L2  =  1227.60  MHz  (24-cm  wavelength)  =  120  x  10.23  MHz. 

The  dual- frequency  operation  is  provided  because  the 
signals  are  delayed  in  passing  through  the  ionosphere  by  an 
amount  inversely  proporational  to  the  frequency  squared.  Mea¬ 
surements  at  two  frequencies  can  be  used  to  compensate  for 
this  ionospheric  delay.  Additional  details  relevant  to  the 
unique  characteristics  and  attributes  of  these  basic  signals 
are  provided  in  Ref.  4,  5,  6,  and  40. 

The  navigation  message  includes  the  following  information 

•  Satellite  status 

•  Satellite  clock  correction  and  ephemeris 
parameters 

•  Ionosphere  propagation  delay  model 
coefficients 

•  Almanac  information  including  the  ephe- 
merides  and  status  of  all  other  satellites. 

This  information  is  employed  in  the  user  equipment  to  make 
"pseudorange”  measurements  to  the  selected  four  satellites. 

The  term  pseudorange  is  used  since  the  measurements  contain  a 
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bias  error  due  to  a  user's  lack  of  a  precise  time  reference 
synchronized  to  GPS  time,  and  is  given  by 

Rt  =  ri  +  c  Atf  +  c(Atu  -  ATi)  (A. 1-1) 

where 

r  v» 

r^  =  true  range  to  i  satellite 
c  =  speed  of  light 

A t^  =  propagation  delays 

At  =  user's  clock  offset  from  GPS  time 
u 

Al\  =  i11*1  satellite's  clock  offset  from  GPS  time. 

GPS  receivers  actually  estimate  both  code  phase  and 
carrier  phase  while  tracking  a  signal  in  their  most  accurate 
mode.  By  differencing  beginning  and  ending  pseudorange  measure 
ments  made  over  a  period  when  carrier  phase  lock  is  achieved, 
an  "integrated  doppler"  measurement  of  "delta  range"  is  possi¬ 
ble.  When  the  doppler  processing  period  is  short  and  the 
delta- range  measurement  is  divided  by  the  period,  the  measure¬ 
ment  is  termed  psuedorange- rate .  Pseudorange-rate  is  the 

time-derivative  of  pseudorange  and  is  a  measure  of  the  rate  at 
which  the  range  from  the  user  to  the  GPS  satellite  is  changing, 
plus  the  frequency  offset  of  the  user  clock.  Just  as  clock 
phase  offsets  or  time  periods  are  measurable  in  units  of  dis¬ 
tance,  so  frequency  offsets  are  measurable  in  units  of  veloc¬ 
ity,  by  multiplying  by  the  speed  of  light. 

Upon  measuring  the  pseudorange  to  four  satellites, 
the  user  equipment  can  compute  the  position  of  these  satel¬ 
lites  using  the  transmitted  ephemeris  data.  Three-dimensional 
user  position  and  precise  time  can  then  be  derived  from  the 
oseudorange  measurements  and  satellite  positions.  This  is  a 
straightforward  matter  involving  the  solution  of  four  simul¬ 
taneous  nonlinear  equations  in  four  unknowns  (see  Ref.  6). 


If  the  errors  in  the  four  simultaneous  pseudorange 
measurements  are  independent  and  of  equal  variance,  then  the 
navigation  errors  resulting  from  the  deterministic  solution 
can  be  related  to  the  pseudorange  errors  by  the  geometric  dilu 
tion  of  precision  (GDOP).  GDOP  refers  to  the  four-dimensional 
root-sum-square  (rss)  combination  of  position  and  time  errors, 

GDOP  =  /aN2  +  a£2  +  aD2  +  (A.  1-2) 

where 

aN’aE’aD  =  1:1113  user  Positi°n  errors  in  a  local- 
A  1  level  North,  East,  Down  coordinate  frame 

at  =  rms  user  time  offset. 

Similar  expressions  exist  for  the  dilution  of  precision  in 
three-dimensional  position  (PDOP),  horizontal  position  (HDOP), 
vertical  position  (VDOP) ,  and  time  (TDOP).  Furthermore,  the 
standard  deviation  of  the  resultant  navigation  error  is  equal 
to  the  appropriate  dilution  of  precision  factor  times  the 
standard  deviation  of  a  single  pseudorange  error  (of) .  Thus, 
the  three-dimensional  rms  position  error  is  given  by 

a  =  a  x  PDOP  (A. 1-3) 

p  r 

The  value  of  PDOP  has  a  geometric  interpretation  as 
well.  Consider  a  hexahedron- shaped  figure  with  the  user  posi¬ 
tion  located  at  the  bottom  indice  and  the  position  of  four 
satellites  located  at  the  other  indices  (see  Fig.  A. 1-3).  The 
PDOP  value  is  inversely  related  to  the  volume  encompassed  by 
this  hexahedron.  Thus,  the  volume  is  maximized  (and  PDOP  is 
minimized)  when  one  satellite  is  at  the  user's  zenith  and  the 
other  three  are  separated  by  120  deg  and  are  as  low  on  the 
horizon  as  the  user's  antenna  elevation  angle  permits  (see 
Ref.  11). 
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Figure  A. 1-3  A  Geometric  Interpretation  of  PDOP 

As  noted  earlier  in  this  section,  during  GGSS  testing 
there  will  be  only  six  satellites  available.  This  situation 
will  lead  to  limited  viewing  times  when  four  satellites  will  be 
visible  simultaneously.  However,  by  equipping  the  user  with  a 
precise  clock,  three-satellite  navigation  can  be  realized  The 
staoility  of  the  clock  is  application-dependent  and  thus,  nec¬ 
essarily  a  function  of  the  maximum  allowable  rms  position  error 
and  clock  update  interval.  Techniques  appropriate  for  updating 
the  clock  include  (Refs.  41  and  42): 

•  Tracking  four  GPS  satellites  and  deriving 
the  nominal  navigation  solution 

•  Tracking  a  single  GPS  satellite  from  a 
known  location  and  extracting  the  appro¬ 
priate  information  from  the  navigation 
message 


•  Conventional  time- transfer  methods. 

As  expected,  the  resultant  TDOP  is  dependent  on  the  accuracy 
of  the  particular  clock  update  technique.  For  example,  for 
the  four-satellite  method, 

TDOP  =  oh/ar  (A. 1-4) 

where  a.  is  the  clock  bias  rms  error  associated  with  the  last 
b 

clock  update.  The  mathematical  basis  for  computing  the  rms 
user  position  error  as  a  function  of  three-satellite  geometry 
and  clock  stability  is  provided  in  Ref.  41. 

The  preceding  discussion  for  the  three-satellite  case 
can  be  naturally  extended  to  apply  to  situations  when  only  two 
satellites  are  visible.  Again,  a  precise  clock  is  necessary, 
but  now  an  altimeter  must  also  be  included  to  provide  the 
vertical  position  information  necessary  for  completing  the 
full  navigation  solution.  Analogous  to  Eq.  A. 1-4, 

VDOP  =  aa/ar  (A. 1-5) 

where  a  is  the  rms  error  of  the  altimeter  measurement, 
a 


A. 2  ORBIT  CHARACTERISTICS 


General  Orbital  Relationship 


Before  discussing  the  orbital  trajectories  of  GPS 
satellites,  it  is  useful  to  review  briefly  the  general  proper¬ 
ties  of  earth  satellite  orbits  (Ref.  43).  Of  the  several 
forces  that  combine  to  shape  such  orbits,  by  far  the  most  im¬ 
portant  is  the  earth's  gravi tational  attraction,  with  lesser 
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contributors  including  atmospheric  drag,  solar  radiation  pres¬ 
sure,  and  gravitational  effects  due  to  the  sun,  the  moon,  and 
the  planets.  The  earth's  gravitational  attraction  can  be 
considered  as  a  sum  of  terms,  with  the  dominant,  central  term 
being  the  simple,  inverse-square,  point-mass  attraction.  If 
this  central  gravitational  term  were  the  only  force  on  a  small, 
relatively  massless  satellite,  then  the  satellite  orbit  would 
be  a  perfect  ellipse  fixed  in  an  earth-centered  inertial  space 
with  one  focus  at  the  center  of  mass  of  the  earth.  Actual 
satellite  orbits  are  nearly,  but  not  exactly,  ellipses,  and 
the  orbits  tend  not  to  remain  fixed  in  inertial  space  over 
long  periods  of  time. 

It  takes,  in  general,  six  parameters  to  describe  the 
state  of  a  point-mass  satellite.  One  such  set  of  parameters 
consists  of  the  vector  components  of  satellite  position  and 
velocity  in  either  an  inertial  or  an  earth- iixed  coordinate 
frame.  Such  a  set  of  parameters  is  very  useful  in  computations 
involving  the  relative  geometry  between  the  satellite  and  other 
points  in  space.  But  it  is  not  a  useful  set  for  conceptualiz¬ 
ing  the  satellite  orbit  or  for  understanding  orbital  propaga¬ 
tion  (evolution),  because  all  six  vector  components  change 
rapidly  in  a  complex  manner.  A  mathematically  equivalent  set 
of  six  parameters  is  more  useful  in  describing  the  satellite 
orbit,  a  set  of  parameters  related  to  the  perfectly  elliptical 
orbit . 


Three  of  the  parameters  describe  the  size  and  shape 
of  the  ellipse  and  the  position  of  the  satellite  in  the  el¬ 
lipse.  These  parameters  are  the  semimajor  axis  length  a,  the 
eccentricity  e,  and  the  (time-varying)  true  anomaly  angle  f, 
all  illustrated  in  Fig.  A. 2-1.  The  remaining  three  angular 
parameters  describe  the  orientation  of  the  ellipse  with  re¬ 
spect  to  the  inertial  coordinate  frame.  They  are  the  orbital 
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inclination  I,  the  right  ascension  of  the  ascending  node  Q, 
and  the  argument  of  perigee  uj  ,  all  illustrated  in  Fig.  A. 2-2. 

It  should  be  noted  that  in  this  appendix  the  right  ascension  of 
the  ascending  node  (the  node  is  the  point  where  the  satellite 
orbit  crosses  the  earth's  equator  going  from  south  to  north) 
is  measured  from  the  x-axis  of  a  basic  inertial  coordinate 
frame  defined  in  Fig.  A. 2-3.  The  orientation  of  this  x-axis 
with  respect  to  the  fixed  stars  is  somewhat  arbitrary,  depending 
on  where  the  Greenwich  meridian  is  at  the  epoch  (tQ)  when  the 
inertial  coordinate  frame  is  conceptually  established  as  coin¬ 
cident  with  earth-fixed  coordinates.  Most  works  (including 
Ref.  44)  measure  the  right  ascension  from  the  vernal  equinox 
(which  is  the  ascending  node  of  the  sun  under  the  fictitious 
concept  that  the  sun  is  orbiting  the  earth). 
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Figure  A. 2-1 
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Figure  A. 2-2 


Orbital  Orientation 


EARTH  s 
SPIN  AXIS 


M 42401 


Figure  A. 2-3 


Inertial  Coordinates 


Also  noted  in  Fig.  A. 2-1  are  coordinate  axes  of  the 
"q"  frame.  The  relationship  between  the  q-frame  and  the  basic 
inertial  frame  is  provided  by  the  three  orientation  parameters , 
and  the  matrix  that  transforms  vectors  expressed  in  the  i- frame 
to  an  expression  in  the  q-frame  is  given  by 


R  .  =  R3(u»)R1(I)R3(n> 


(A. 2-1) 
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(A. 2-2) 


where  the  C  and  S  designations  represent  the  cosine  and  sine 
trigonometric  functions,  respectively. 

The  motion  of  the  satellite  in  its  elliptical  orbit 
will  not  be  derived  (see  Ref.  44),  but  key  results  will  be 
summarized.  The  position  vector  from  the  center  of  mass  of 
the  earth  (the  focus  of  the  elliptical  orbit)  to  the  satellite 
is  given  by 


"  r  Cr 
S  f 

( r  )  =  r  S  f 

s  q  s  f 


(A. 2-3) 


where  rg  is  the  magnitude  of  the  radius  vector.  In  terms  of 
the  elliptical  parameters 


_  a  ( 1  -  e  ) 


r 


(A. 2-4) 


Of  the  six  parameters  describing  the  satellite  in  its 
idealized  elliptical  orbit,  five  are  constants,  and  only  one, 
the  true  anomaly  (f)  changes  with  time.  The  way  in  which  the 
true  anomaly  changes  with  time  is  somewhat  complex,  and  is 
best  described  by  defining  two  new  terms:  the  eccentric  anomaly 
E,  and  the  mean  anomaly  M.  All  of  the  three  anomalies  are 
defined  to  be  zero  when  the  satellite  is  at  perigee. 

The  mean  anomaly  changes  at  a  constant  rate 

M  =  n(t-t  )  (A. 2-5) 

where  t  is  the  time  of  perigee,  and  n  is  the  mean  motion  to 

r 

be  discussed  later.  The  eccentric  anomaly  is  related  to  the 
mean  anomaly  via  Kepler's  Equation 

E  -  eSg  =  M  (A. 2-6) 

The  relationship  between  the  eccentric  and  the  true  anomalies 
is  best  provided  by  the  equations 

rsCf  =  a^CE’e>  (A-2-7) 

rgSf  =  a  Jl-e2  (A. 2-8) 


A  geometric  justification  for  the  definition  of  the  eccentric 
anomaly  is  provided  in  Ref.  43.  An  important  point  to  note 
here  is  that  if  the  satellite  orbit  is  circular  then  the  eccen¬ 
tricity  e  is  zero,  and  the  true  anomaly,  the  eccentric  anomaly, 
and  the  mean  anomaly  are  all  identical. 


The  rate  of  the  mean  anomaly,  the  mean  motion  (n),  is 
related  to  the  earth's  gravitational  constant  in  the  assumed 
central  force  term  (Mg)  and  to  the  orbit  semimajor  axis  (a)  as 
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The  mean  motion  shows  up  in  the  expression  for  the  satellite 
(inertial)  velocity  vector  in  the  q-frame 
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Anticipating  the  fact  that  GPS  satellite  orbits  are 
nominally  circular,  the  foregoing  expressions  can  be  simplied 
in  two  ways.  The  first  way  is  by  setting  the  orbital  eccen¬ 
tricity  to  zero.  The  second  simplification  comes  from  the  fact 
that  there  is  no  perigee  in  a  circular  orbit,  so  the  arbitrary 
argument  of  perigee  can  also  be  set  to  zero,  and  the  various 
anomalies  (now  all  identical)  can  be  measured  from  the  ascend¬ 
ing  node.  For  a  reminder,  denote  the  anomalies  all  by  the 
symbol  A.  The  desired  expressions  are  the  satellite  position 
and  velocity  components  in  the  inertial  coordinate  frame  in 
terms  of  the  elliptic  orbit  parameters.  They  are: 
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(A. 2-11) 


(A. 2-12) 


(A. 2-13) 


(A. 2-14) 
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The  equivalent  expressions  in  earth-fixed  coordinates  are  iden¬ 
tical  except  that  the  right  ascension  Q  must  be  replaced  by 
Q ' ,  where 

Q'  =  Q  -  eg  (A. 2-15) 

9o  being  the  earth's  sidereal  angle  defined  by 

9e  =  we(l:  '•  to)  (A. 2-16) 

where  uie  denotes  the  angular  velocity  of  the  rotating  earth 
with  respect  to  inertial  space. 

The  next  step  in  this  argument  is  to  connect  idealized 
elliptical  orbits  with  real,  nearly  elliptical  orbits.  The 
connection  is  provided,  on  an  instantaneous  basis,  by  the 
relationship  between  the  elliptical  orbit  parameters  and  the 
satellite  position  and  velocity  components  in  inertial  space 
just  derived.  At  any  instant  in  time  the  satellite  has  some 
definite  position  and  velocity  components  in  the  basic  inertial 
frame.  There  is,  at  the  same  instant,  some  set  of  elliptical 
orbit  parameters  that  correspond  to  the  same  position  and  veloc 
ity  components  via  the  relationships  above.  This  elliptical 
orbit  is  termed  an  instantaneous  osculating  orbit.  Thus,  any 
real  orbit  (of  the  nearly  elliptical  category  --  as  opposed  to 
the  nearly  hyperbolic)  can  be  defined  in  terms  of  a  set  of 
elliptical  orbit  parameters.  In  an  idealized  elliptical  orbit, 
five  of  the  six  orbit  parameters  are  constants,  while  one 
changes  with  time.  In  a  real  orbit  described  by  elliptical 
orbit  parameters,  all  six  parameters  are  subject  to  change 
with  time.  The  advantage  of  using  the  elliptical  parameters 
is  that  five  of  the  parameters  should  change  very  slowly. 

The  chief  disturbance  of  the  elliptical  orbits  of  most 
high-altitude  satellites  (those  of  reasonably  high  density)  is 


che  terra  in  the  earth's  gravi tational  field  resulting  from  the 
ellipsoidal  figure  of  the  earth.  The  coefficient  of  this  term, 
often  noted  as  C 2q  (or  is  over  100  times  as  great  as  the 

coefficient  of  any  other  (except  the  central)  term  in  the 
spherical  harmonic  expansion  of  the  earth's  gravitational 
potential.  The  effect  of  the  earth's  equatorial  bulge  on 
satellite  orbits  can  be  computed  approximately  by  orbital 
perturbation  theory.  Ignoring  small  cyclic  effects  during  a 
single  orbit  and  concentrating  on  long-term  trends,  the  results 
are  that  the  bulge  does  not  alter  the  orbit  semimajor  axis, 
the  eccentricity,  or  the  inclination  angle,  but  does  cause 
trends  in  the  argument  of  perigee  and  the  right  ascension, 
and  alters  the  value  of  the  mean  anomaly  rate  slightly.  Direc¬ 
tly  from  Ref.  44,  simplified  for  a  circular  orbit,  the  per¬ 
turbed  results  are 
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It  might  seem  strange  to  express  a  rate  of  change  of 
the  argument  of  perigee  for  a  circular  orbit  that  has  no  peri¬ 
gee  and  for  which  it  was  previously  assumed  that  the  argument 
of  perigee  was  zero.  A  better  way  to  think  about  the  previous 
manipulations  in  this  perturbed  case  is  that  the  argument  of 
perigee  is  initial ized  at  zero,  and  the  "common"  anomaly  (A) 
is  actually  the  sum  of  the  argument  of  perigee  and  the  mean 
anomaly  (or  the  eccentric  anomaly,  or  the  true  anomaly).  Now 


A  =  M  +  ui 


(A. 2-20) 
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The  values  of  the  two  gravitational  constants  involved 
in  the  above  relationships,  the  gravitational  constant  pe  and 
the  second  degree  zonal  coefficient  C2 q,  are  defining  param¬ 
eters  of  the  WGS  72  ellipsoid,  and  are 


Me  =  3.986005  x  1014  m3/sec2 
C2Q  =  -4.841605  x  10‘4 


The  unnormalized  coefficient  used  in  Ref.  44  and  in  the  above 
equations  is  related  to  the  normalized  version  by 

C2Q  =  J5  C2Q  (A. 2-22) 

so  that 

C2Q  =  -1.0826158  x  10"3 


Because  this  value  is  negative,  it  is  apparent  from  Eq.  A. 2-18 
that  all  orbits  with  inclination  angles  less  than  90  deg  will 
have  a  decreasing  right  ascension.  That  is,  their  orbit  planes 
will  rotate  slowly  in  a  direction  opposite  to  the  earth's 
rotation  rate.  To  proceed  further,  it  is  necessary  to  know 
more  about  the  nominal  orbit  parameters,  so  this  discussion 
will  be  suspended  until  some  details  of  the  GPS  satellite 
orbits  are  introduced. 


A . 2 . 2  Common  Features  of  GPS  Satellite  Orbits 

In  the  Global  Positioning  System  (GPS),  essentially 
line-of-sight  (L-band)  navigation  signals  are  broadcast  by  a 
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number  of  earth  orbiting  satellites.  The  set  of  satellites  is 
termed  the  GPS  constellation.  Specific  characteristics  of  the 
current  test  and  planned  operational  GPS  constellations  are 
identified  in  Section  A.l. 

All  of  the  GPS  satellites  will  be  in  nominally  circu¬ 
lar  "12-hour"  orbits.  The  common  orbital  radius  of  the  satel¬ 
lites  is  chosen,  as  will  be  shown,  to  yield  an  orbital  period 
of  very  nearly  one-half  of  a  sidereal  day.  A  sidereal  day  is 
the  time  required  for  the  earth  to  complete  one  full  rotation 
with  respect  to  the  fixed  stars,  and  since  the  satellites 
complete  two  orbital  revolutions  in  the  same  time,  the  posi¬ 
tions  of  the  satellites  relative  to  the  earth  are  periodic, 
with  a  period  of  very  nearly  one  sidereal  day  (one  ordinary 
day  less  about  four  minutes).  The  "very  nearly"  qualifier  is 


added  because,  as  was  noted  in  the  previou  section,  the  satel 
lite  orbits  are  not  quite  fixed  with  respect  to  the  stars,  but 
the  orbital  planes  rotate  slowly  due  primarily  to  the  oblate¬ 
ness  of  the  earth. 

The  GPS  satellites  will  all  have  the  same  nominal 
orbital  inclination  of  55  deg  (it  is  63  deg  in  the  current 
test  constellation).  Using  the  approximate  perturbation  re¬ 
sults  of  the  previous  section,  it  is  possible  to  compute  the 
nominal  radius  and  period  of  the  GPS  satellite  orbits.  The 
key  equation  can  be  deduced  from  the  fact  that  the  satellite 
ground  tracks  must  be  periodic.  Thus,  Eqs.  A. 2-12  and  A. 2-14 
(with  Q'  replacing  Q  so  that  they  refer  to  the  earth-fixed 
frame  rather  than  the  inertial  frame)  must  be  periodic.  This 
can  only  happen  if  the  anomaly  A  and  the  earth  relative  "right 
ascension"  Q'  share  a  common  period.  For  the  "12-hcur"  orbits 
the  anomaly  rate  must  be  twice  the  effective  rate  of  the  earth 
rotation  with  respect  to  the  GPS  satellite  orbit  planes.  In 
symbols , 


(A. 2-23) 


A  =  -2  Q'  =  2(6e  -  Q) 

Upon  substituting  Eq.  A. 2-21  for  A  and  Eq .  A. 2-18  for 
Q,  recalling  that  the  rate  of  change  of  the  sidereal  angle  0g 
is  just  the  sidereal  earth  rate  wg  ,  using  Eq .  A. 2-9  for  the 
mean  motion  n,  and  rearranging  yields 


(A. 2-24) 

Everything  in  this  relationship  is  known  except  the  GPS  satel¬ 
lite  radius  r  .  Numerically  solving  for  the  ratio  <ra/aj 

S  S  0 

yields 


(rs/ae)  =  4.1642875 

or 

r  =  26560387  m 
s 

Having  the  GPS  satellite  allows  some  other  relevant 
Quantities  to  be  computed.  The  satellite  orbital  rate  is 

A  =  1.4585797  x  10"4  rad/sec 
the  orbital  period  is 

TGps  =  2n/k  =  llh  57m57 . 42s  (A. 2-25) 

-the  ground  track  period  is 

2Tcps  =  23h  55m  54.84s 
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the  right  ascension  rate  is 


Q  =  -7.834  x  10  9  rad/sec 
Q  =  -Id. 17  deg/yr 
the  mean  motion  is 

n  =  1.4585366  x  10  ^  rad/sec 
and  the  satellite  orbital  velocity  is 
vg  =  3873.93  m/sec 

An  extensive  computer  simulation,  accounting  for 
numerous  disturbing  effects,  beyond  the  major  one  of  the  earth's 
obiateness,  would  result  in  values  for  the  key  parameters  that 
differ  very  little  from  those  presented  above. 

A. 2. 3  Simulating  the  GPS  Satellite  Constellation 

The  longitude  and  time  connections  between  the  opera¬ 
tional  GPS  satellite  constellation  and  earth-fixed  coordinates 
are  not  yet  finalized.  In  order  to  parameterize  these  connec¬ 
tions  and  understand  them  more  fully,  it  is  useful  to  introduce 
a  new  parameter:  the  longitude  of  the  ascending  node.  Up  co 
this  point,  the  location  of  the  ascending  node  of  a  satellite 
orbit  has  been  described  by  its  right  ascension  angle  measured 
from  the  x^  axis  of  the  basic  inertial  frame  (the  location  of 
the  Greenwich  meridian  at  some  arbitrary  time  tQ).  The  right 
ascension  parameter  is  most  useful  when  dealing  with  idealized 
elliptic  orbits  where  the  right  ascension  does  not  vary  with 
Lime.  For  GPS  satellite  orbits,  whose  most  constant  aspect  is 
their  repeating  ground  tracks,  the  longitude  of  the  ascending 
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node  is  more  useful.  It  is  defined  to  be  the  satellite  longi¬ 
tude,  measured  from  the  Greenwich  meridian,  at  the  instant 
when  the  satellite  is  aj:  its  ascending  node.  Since  the  satel¬ 
lite  is  at  its  ascending  node  once  each  12  hours  or  so,  there 
is  a  sequence  of  longitudes  of  the  ascending  nodes.  But  for 
GPS  satellites,  this  sequence  is  just  two  values,  differing  by 
180  deg,  repeated  alternately.  One  of  these  two  values,  the 
one  greater  than  or  equal  to  zero  and  less  than  180  deg,  can 
arbitrarily  be  selected  as  the  primary  longitude  of  the  ascend 
ing  node.  The  longitude  of  the  ascending  node  for  the  #1  GPS 
Satellite  is  used  to  parameterize  the  longitude  connection 
between  the  satellite  constellation  and  the  earth,  and  is 
denoted  here  by  \Gpg. 

The  time  connection  parameter  is  the  time  of  the 
ascending  node  of  the  #1  GPS  Satellite,  denoted  by  tGps .  Of 
course,  there  is  a  sequence  of  such  times,  separated  by  ap¬ 
proximately  24  hours,  since  it  is  the  times  of  the  primary 
nodes  that  have  longitude  AGps  that  are  being  considered.  Any 
member  of  the  sequence  will  serve  to  establish  the  time  connec 
tion.  It  is  often  easiest  to  use  that  member  closest  to  the 
"time  zero"  that  establishes  the  basic  inertial  reference 
system  with  respect  to  earth-fixed  coordinates. 

In  terms  of  these  new  parameters,  the  right  ascension 
of  the  ascending  node  of  the  #1  GPS  Satellite,  at  the  time 
tGps,  is  given  by 


°GPS  ~  XGPS  +  We(tGPS  “  to) 


(A. 2-26) 


At  any  time  t,  the  right  ascension  of  the  j  GPS  satellite  is 
given  by 


0  •  ( t )  =  Q„np 


+  Q  „  .  +  fl(t-t„„*.) 


(A. 2-27) 
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where  is  the  relative  right  ascension  of  the  j  satellite 

(tabulated  in  Table  A. 1-2  or  otherwise  known  from  the  constel¬ 
lation),  and  Q  is  the  constant,  very  small,  right  ascension 
rate  given  following  Eq.  A. 2-25. 

The  anomaly  of  the  j t^1  GPS  satellite  at  any  time  t  is 
given  by 


Aj  ^ Arj  +  A^  t_tGPS  ^ 


(A. 2-28) 


where  A  .  is  the  relative  anomaly  (tabulated  in  Table  A. 1-2  or 
rj  3 

otherwise  known  from  the  constellation),  and  A  is  the  not-so- 
small,  constant  value  given  just  prior  to  Eq .  A. 2-25. 
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Table  A. 2-1  summarizes  the  data  and  calculations  nec¬ 
essary  to  simulate  a  GPS  constellation  in  terms  of  the  posi¬ 
tions  and  velocities  of  its  satellites  in  the  basic  inertial 
frame.  For  short  time  period  simulations,  the  propagation 
equation  for  the  satellite  right  ascensions  can  be  used  to 
initialize  the  right  ascension  values  which  can  subsequently 
be  treated  as  constants,  avoiding  the  reevaluation  of  numerous 
trigonometric  functions  as  time  passes.  Further,  for  a  fixed 
time-step  simulation,  the  sine  and  cosine  addition  formulas 
can  be  used  to  propagate  the  sines  and  cosines  of  the  satellite 
anomalies,  again  avoiding  the  trigonometric  function  computa¬ 
tions.  That  is,  over  a  fixed  time  step  At,  each  satellite 
anomaly  advances  by  an  amount 


AA  =  A  At 


(A. 2-29) 
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and  the  sine  and  cosine  functions  advance  as 


S(A.+AA)  =  SA.CAA  +  CA.SAA 
J  J  J 


(A. 2-30) 


C ( A . +AA)  =  CA.CAA  "  SA.SAA 
J  J  J 


(A. 2-31) 


Thus,  the  sine  and  cosine  of  the  fixed  anomaly  step  can  be 
computed  but  once,  as  can  the  sines  and  cosines  of  the  initial 
anomalies  for  all  the  satellites.  Subsequently,  four  multi¬ 
plies  and  two  adds  serve  to  propagate  the  sine  and  cosine  of 
each  satellite  anomaly.  Another  possible  simplification  of 
the  simulation  computations  comes  from  the  fact  that  the  satel 
lite  anomalies  of  the  various  satellites  are  sometimes  equal 
or  separated  by  integral  quadrants,  making  for  especially 
simple  relationships  among  the  sines  and  cosines. 


Besides  computing  the  position  and  velocity  vectors 
of  the  GPS  satellites,  it  is  sometimes  desired  to  compute  the 
acceleration  and  jerk  (acceleration  rate)  vectors.  Because  of 
the  circular  orbits,  the  acceleration  and  jerk  are  related  to 
the  position  and  velocity  by 


.  2 

a  =  -A*  r 

(A. 2-32) 

— s  — s 

•  2 

j  =  -A  v 
“-S  -s 

(A. 2-33) 

A. 3  ANALYSIS  RESULTS 

Satellite  visibility  is  a  key  issue  which  must  be 
investigated  when  analyzing  the  ability  of  GPS  to  satisfy  GGSS 
positioning  and  navigation  requirements.  In  this  section, 
various  visibility  factors  are  considered  and  analysis  results 
are  presented  for  the  complete  18-satellite  constellation. 


The  orbit  relationships  of  the  operational  GPS  constel¬ 
lation  are  presented  in  Section  A.l.  Current  plans  call  for 
the  GPS  antenna  to  be  mounted  on  the  top  of  the  aircraft  fusel¬ 
age.  In  such  an  arrangement,  the  earth  is  not  a  factor  in 
signal  blockage.  In  addition,  while  the  gravity  gradient  data 
are  collected  during  an  individual  survey  sortie  (see  Fig.  A-3-1 
for  a  typical  pattern),  the  antenna  field  of  view  is  not  limited 
As  the  plane  turns  and  banks  to  begin  another  data-gathering 
leg  of  the  sortie,  the  wing  can  block  the  signal  temporarily. 
However,  this  condition  should  not  matter  since  the  position¬ 
ing  data  are  not  needed  in  this  area  beyond  the  survey  region 
and  the  GGSS  navigation  system  dynamics  are  being  designed  to 
recover  quickly  from  a  situation  such  as  this. 

Using  the  relative  Block  II  orbit  parameters  listed 
in  Table  A. 1-2,  a  GGSS  user's  latitude,  longitude,  and  nominal 
600  m  altitude,  an  appropriate  user  elevation  mask  angle  of 

5km 
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Figure  A. 3-1  A  Proposed  GGSS  Survey  Sortie 
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Typical  Visibility  Pattern  for 
18-Satellite  Constellation 
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However,  to  completely  characterize  the  effectiveness 
of  using  GPS  as  a  GGSS-aiding  sensor,  the  geometry  associated 
with  the  satellite  visibility  pattern  must  also  be  determined. 
Figure  A. 3-3  is  a  plot  of  the  PDOP  which  is  associated  with 
the  pattern  depicted  in  Fig.  A. 3-2.  Since  there  were  several 
instances  when  more  than  four  satellites  are  visible,  the  PDOP 
which  reflects  the  selection  of  the  four  "best"  satellites 
(i.e.,  those  which  result  in  the  minimum  value  of  PDOP)  is 
plotted.  The  statistics  (mean  and  standard  deviation)  of  this 
curve  are  also  indicated  on  the  figure.  Recall  that  the  PDOP 
is  directly  proportional  to  the  three-dimentional  rms  position 


error ,  i . e . , 


a  =  PDOP  x  a 
P  r 


(A. 3-1) 


MEAN  -  2  80 


STANOARO 
DEVIATION  =  0  70 


where  is  the  rms  user-to-satellite  range  error  which  is 

commonly  assumed  to  be  7  m.  This  value  is  consistent  with 
test  results  to  date  (Refs.  38  and  39).  Thus,  for  the  case 
depicted  in  Fig.  2.2-3,  ~  20  m  rms  which  is  well  within 

the  100  m  rms  GGSS  positioning  accuracy  goal  identified  in 
Section  2.1. 

As  pointed  out  in  Ref.  39,  a  user  will  experience  short 
periods  of  poor  accuracy  even  with  the  complete  18-satellite 
constellation  operational.  These  "outages"  are  due  solely  to 
poor  geometry  since  at  least  four  satellites  will  always  be 
in  view.  Only  one  set  of  four  widely  separated,  predictable 
outages  exists  at  any  time  and  the  time  duration  of  an  outage 
at  any  of  the  particular  locations  is  quite  short  (e.g.,  5  to 
20  min).  The  outages  repeat  twice  a  day  at  these  locations, 
disappear,  and  then  appear  40  min  later  at  four  other  locations. 

A  composite  of  all  outages  which  will  occur  during  a 
given  day  is  presented  in  Fig.  A. 3-4.  For  an  arbitrary  start¬ 
ing  time,  the  set  of  four  labeled  1  occurs  first,  followed 
40  min  later  by  the  set  designated  2,  etc.  Within  the  circled 
outage  regions,  the  PDOP  will  be  very  large  (on  the  order  of 
thousands)  and  result  in  unacceptable  four-satellite  navigation 
accuracy  levels.  A  plot  of  the  PDOP  for  a  user  at  one  of  these 
outage  locations  is  presented  in  Fig.  A. 3-5.  As  is  apparent 
from  this  figure  and  the  included  PDOP  statistics,  since  the 
outage  intervals  are  predictable  and  short,  they  are  inconse¬ 
quential  to  the  currently  configured  GGSS. 
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Figure  A. 3-4  Composite  Outages  for  the  18-Satellite 

Constellation  (Extracted  from  Ref.  39) 
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APPENDIX  B 


DERIVATION  OF  DOPPLER- SHIFT  OBSERVABLE  AND  ASSOCIATED 
VEHICLE  ERROR  DYNAMICS  FOR  APPLICATIONS  OF 
MOVING-RECEIVER  RADIO  INTERFEROMETRY 


B.l  INTRODUCTION 

This  appendix  details  the  derivations  of:  1)  the 

doppler-shift  measurement  equation  (Eq.  3.2-5)  and  2)  an  ex¬ 
pression  for  the  vehicle  velocity  error  (Eq.  3.4-2).  Both 
equations  depend  on  the  basic  doppler-shift  definition 


f ' 


f 


1  +  nvcosQ/c 
1  -  nvcosQ/c 


(B.l-1) 


where 

f'  =  doppler-shifted  frequency  observed  at  the 
moving -receiver 

f  =  frequency  observed  at  the  fixed-site 
receiver 

n  =  index  of  refraction 
c  =  free-space  light  speed 
v  =  moving- receiver  velocity 
0  =  satellite  signal  incidence  angle. 


B . 2  DERIVATION  OF  DOPPLER-SHIFT  OBSERVABLE 

Equation  B.l-1  can  be  easily  expanded  using  a  Mac- 


laurin  series: 


(B.  2-1  ) 


z(x)  =  z(x  =  0)  +  X 


dz 

dx 


_  X  d  Z 

2:  .  2 

A  QX  I  .fl 

x=0  x=0 


Expanding  Eq .  B.l-1  in  terms  of  y  =  v/c  yields 


f ' (y  =  0)  =  f 


(B.2-2 


df ' 
dy 


_ fncosQ _ 

( l-ycos8  )^//2  ( 1 +y  c  o  s  0  )^//^ 


(B.2-3) 


d  f '  _  fn  cos  0(l+2ycos0) 

— T~  ~  - - - T77 

dy  ( 1  -ycose  )  z  ( l-*-ycos0  ) 


(B.2-4) 


Substitution  of  Eqs .  B.2-2,  B.2-3,  and  B.2-4  into  B.2-1  leads 
to 

f'  =  f [ 1  +  nvcos6/c  +  1/2 (nvcos0/c ) ^ ]  +  0(10  (B.2-5) 


Finally,  the  difference  in  the  doppler- shi f t ,  Af,  as  recorded 
at  moving-  and  fixed-site  receivers,  can  be  written  as 


Af  =  f’ 


_  fnvcosO  [ i  ,  nvcos0 
c  L1  —2C— 


(B.2-6) 


Thus,  an  expression  for  the  doppler- shi ft  difference,  which  is 
one  of  the  observables  associated  with  moving- rece iver  radio 
interferometry,  has  been  derived.  Note  that  vehicle  velocity 
expl ic i t ly  appears  in  this  expression. 


B . 3  DERIVATION  OF  VEHICLE  ERROR  DYNAMICS 

Equation  B.l-1  can  be  inverted  to  obtain  the  followm 
relationship  tor  vehicle  velocity  in  terms  of  frequency: 


Following  a  first-order  perturbation  of  Eq .  B.3-1,  the  error 
in  velocity  can  be  written  as 


5v 


=  -  <H> 


6n  +  vtan6  60  + 


4ff 


ncos0 ( f 


7T 


f2) 


6  f 


(B. 3-2) 


Note  that  the  last  term  in  Eq .  B.3-2  can  also  be  written  in 
terms  of  velocity;  hence 


6v  =  6p  -  - 


<S> 


6n  +  vtan0  60  + 


2  2 

v ( c  - ( nvcosQ ) 

f  ( cncos0 ) 


5  f 


(B.3-3) 


2  2 

Letting  4  =  1 c  "  (nvcos0)  ]/(cncos0),  Eq .  B.3-3  can  be  re¬ 
written  as 

6p  =  -  (-)  6n  +  vtan0  68  +  (4^)  6f  (B.3-4) 

n  t 


Equation  B.3-4  provides  the  relationship  between  the  time  rate 
of  change  of  position  error  (6p),  vehicle  velocity,  and  meas¬ 
ured  values  of  re f rac tivi ty ,  incidence  angle,  and  signal 
frequency . 


APPENDIX  C 

ATMOSPHERIC  PROPAGATION  DELAY 


C . 1  BACKGROUND 

Passage  of  electromagnetic  waves  through  the  atmos¬ 
phere  is  described  by  the  branch  of  wave  mechanics  known  as 
geometrical  optics.  In  this  field  of  science,  electromagnetic 
energy  is  transported  along  light  rays  and  the  velocity  propa¬ 
gation  of  these  rays  is  described  by 

v  =  c/n  (C.l-1) 

where  n  is  the  index  of  refraction  and  c  is  the  free-space 

O 

light  speed  (2.9978  x  10  m/sec).  Depending  on  the  dielectric 
properties  of  the  medium,  the  refractive  index  is  greater  than 
or  less  than  unity. 

The  dielectric  properties  of  the  atmosphere  produce 
refractive  indices  which  are  very  region-dependent.  The  lower 
40  km  (called  the  troposphere)  consists  of  isotropic,  neutral 
gases  which  yield  a  refractive  index  greater  than  unity  and 
which  is  independent  of  wave  polarization  or  propagation  direc¬ 
tion.  The  upper  atmosphere  (collectively  termed  the  ionosphere) 
comprises  anisotropic,  charged  gases  which  produce  two  refrac¬ 
tive  indices  (both  less  than  unity)  and  two  characteristic 
polarizations.  The  values  of  the  indices  and  polarizations 
depend  on  the  wave  propagation  direction. 

A  change  in  the  atmospheric  refractive  index  of  only 
-  few  parts  per  million  can  have  a  significant  effect  on  radio¬ 
waves.  Since  refractive  values  are  very  near  unity  (typically 


n  =  1.00035)  meteorologists  have  defined  the  term  re f rac t ivi ty  , 
N ,  to  be 


N  =  (n-1)  x  106  (C . 1-2 ) 

The  means  by  which  atmospheric  properties  are  measured 
to  yield  refractivity  is  the  subject  of  this  Appendix.  Sec¬ 
tion  C.2  describes  the  methods  employed  in  recording  tropospheric 
refractivity;  Section  C.3  details  the  physics  of  ionospheric 
signal  delay;  and  Section  C.4  describes  methods  to  correct  for 
ionospheric  refraction. 


C.2  TROPOSPHERIC  REFRACTIVITY 

It  is  in  the  troposphere  that  changes  of  temperature, 
humidity,  and  pressure,  as  well  as  clouds  and  rain,  influence 
the  way  in  which  radio  waves  propagate.  Spot  measurements  of 
the  tropospheric  refractive  index  are  made  using  a  variety  of 
application-dependent  techniques . 

For  airborne  GGSS  survey  monitoring,  ruggedness, 
repeatability  without  frequent  calibration,  and  a  wide  range 
of  sensor  response  are  more  important  than  high  accuracy.  If 
a  system  is  needed  to  measure  rapid  index  fluctuations  (at 
a  rate  of  several  tens  of  Hertz),  only  instruments  which  record 
the  refractive  index  directly  have  an  adequate  speed  of  re¬ 
sponse;  such  instruments  are  referred  to  as  "ref ractometers . " 

If  a  lower  sampling  rate  can  be  permitted,  or  a  lightweight 
system  is  essential,  the  refractive  index  can  be  determined 
indirectly  through  measurements  of  atmospheric  properties. 

Microwave  re  f ractometers  are  the  most  common  means  of 
measuring  the  re  tractive  index  directly. 


The  principle  of 


operation  is  to  measure  the  change  in  the  resonant  frequency, 
6f,  of  a  cavity  with  partly  open  ends.  This  change  occurs 
due  to  a  variation  in  the  refractive  index  of  the  air  passing 
through  the  cavity.  Then,  for  a  fixed  cavity,  the  index  vari¬ 
ation,  6n,  is 

6n  =  -6  f/f  (C.2-1) 

Refractometers  have  been  used  successfully  since  1950  in  air¬ 
crafts,  on  tall  masts,  and  suspended  from  tethered  balloons. 
With  solid-state  circuitry  and  solid-state  microwave  sources, 
the  weight  of  the  sensors  has  been  substantially  reduced  and 
the  short-term  errors  have  been  lowered  to  typically  less  than 
0.01  N-unit.  Stability  of  one  part  in  10^  per  month  has  been 
reported  (Ref.  45). 

Measurement  of  temperature,  water  vapor  content,  and 
air  pressure  can  also  be  used  to  determine  refractivity  in¬ 
directly.  The  empirical  formula  relating  these  quantities  to 
the  refractive  index  is  (Ref.  20) 

N  =  77.6  P/T  +  3.73  *  105  e/T2  (C.2-2) 

where  P  is  the  total  air  pressure  (in  mbar) ,  e  is  the  partial 
pressure  of  water  vapor  (in  mbar),  and  T  is  the  temperature 
(in  Kelvin).  This  equation  has  an  inherent  accuracy  of  ±0.5 
percent  for  atmospheric  pressures  between  200  and  1100  mbar, 
air  temperatures  between  240  and  310  K,  water  vapor  pressures 
less  than  30  mbar,  and  radio  frequencies  less  than  30  GHz. 
For  GPS  applications  (with  a  maximum  frequency  of  1.5754  GHz), 
typical  atmospheric  conditions  near  ground  are  P  =  1000  mbar, 
T  =  288  K,  and  e  =  11.9  mbar  (which  implies  70%  humidity). 

Conventional  radiosondes ,  as  used  for  weather  fore¬ 
casting,  have  been  developed  as  moderately  robust,  low-cost, 
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throw-away  devices  to  measure  pressure,  temperature,  and  humidity 
at  heights  up  to  40  km.  However,  the  sensitivities  and  rapid¬ 
ity  of  response  are  often  inadequate  for  radio-meteorological 
studies.  This  has  prompted  the  development  of  specialized 
radiosondes  which  are  suspended  below  tethered  or  free  bal¬ 
loons.  Many  of  these  sondes  have  employed  direct  measurement 
of  dewpoint  temperatures  to  achieve  rms  refractivity  errors  as 
low  as  1  N-unit,  with  1.3-second  time  constants  (Ref.  20). 

An  important  point  to  remember  about  radiosonde  meas¬ 
urements  of  the  refractive  index  is  the  local  nature  of  the 
recordings;  both  types  of  sensors  record  the  atmosphere  in  the 
vicinity  of  the  monitor.  However,  refractivity  varies  with 
altitude  as  well  as  with  horizontal  distance  separation.  Until 
recently,  this  contribution  to  tropospheric  signal  delay  was 
compensated  only  by  approximate  modeling.  Water-vapor  radio¬ 
meters  now  exist  which  can  measure  the  delay  time  directly. 

Dual- frequency  microwave  radiometers  measure  the  ther¬ 
mal  atmospheric  emission  at  22.2  GHz  (for  water  vapor  resonance) 
and  31.4  GHz  (for  liquid  water  contribution)  to  determine  the 
1 ine-of-sight  signal  path  delay.  These  instruments  have  been 
used  successfully  in  the  Crustal  Dynamics  Program  (Ref.  26) 
and  in  NGS  tests  at  Boulder,  Colorado  (Ref.  27).  While  current 
radiometers  are  not  portable,  ongoing  research  at  JPL  involves 
developing  a  lightweight,  easily  movable  version.  A  prototype 
version  is  anticipated  soon. 

C . 3  IONOSPHERIC  REFRACTIVITY 

The  refractive  index  of  the  upper  atmosphere  is  a 
function  of  radio  frequency  and  free-elec tron  density.  Using 
magne to- ionic  theory  (Ref.  46)  a  relationship  can  be  derived 
between  n  and  the  radio  frequency  f: 


(C. 3-1) 


nA>r  =  U  -  fp2/(f2  ±  f fg ! cose ! ) ] 1/2 


where 


P 

A 


N 


"g 


B 


plasma  frequency;  f  =  ANg 
80.6  m3/sec2 

3 

number  density  of  free  electrons  (m  ) 

frequency  of  gyration  (gyro  frequency)  of  free 

electrons;  f  =  eB/(2nm) 

g 

2 

geomagnetic  induction  (v-sec/m  ) 

angle  between  the  phase  propagation  and  geomagnetic 
induction  vectors. 


The  positive  sign  corresponds  to  the  left-handed  ( 1  - subsc r ip  ted ) 
circular  component  of  polarization  in  the  northern  geomagnetic 
hemisphere.  Expanding  the  square  root  in  Eq .  C.3-1  leads  to 

n2>r  =  1  -  fp/ ( 2f 2 )  ±  f2fg|cos0|/f3  -  f4/(8f4) 

-  f2 f2 (sin2  +  2cos28 )/(4f4)  ±  •••  (C.3-2) 

P  g 

which  is  complete  to  order  (1/f4). 

Equation  C.3-2  can  be  separated  into  angular  and  non- 
angular  dependent  terms.  The  two  angular-dependent  terms 
account  for  the  effect  of  wave  propagation  in  the  geomagnetic 
field  and  are  functions  of  electron  content  and  signal  polari¬ 
zation.  (Circular  polarization  yields  the  extreme  positive  or 

3 

negative  values  of  the  1/f  term. ) 

J 

4 
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Plasma  modeling  of  the  ionosphere  is  required  to 
evaluate  accurately  Eq.  C.3-2.  This  modeling  requires  defini¬ 
tion  of  the  "effective  electron  thickness,"  t.  In  terms  of 


the  electron  density,  N  , 

m 


t  =  /  Ne  dh/Nm  =  I/Nm 


(C. 3-3) 


where  hg  is  the  altitude  of  the  satellite,  I  is  the  total 


electron  content,  and  Ng  is  the  number  density  of  free  electrons 


This  last  parameter  is  a  function  of  geographic  and  geomagnetic 


coordinates,  solar  activity,  season,  time  of  day,  and  geomag¬ 


netic  activity;  Ng  varies  from  1  x  lO^^m”^  to  2  x  10^m-2.  The 


uniformity  or  "shape"  of  the  electron  content  can  be  modeled  by 


1  5  Ne  2 

y  =  t  l  <ir> dh 

u  m 


(C. 3-4) 


In  general,  it  is  not  possible  to  gain  accurate  information  on 
the  height  distribution  of  free  electrons  and  the  "worst-case" 
scenario  implies  a  shape  factor  of  y  =  1,  i.e.,  a  uniform 
electron  density  throughout  the  ionosphere  with  a  thickness  of 
i  =  200  km  (Ref .  47 ) . 


The  delay  of  electromagnetic  signals  propagating 
through  the  ionosphere  is  directly  related  to  the  medium's 
refractive  index.  For  radio  signals  of  f  =  2  GHz  propagating 
along  the  direction  of  geomagnetic  induction  (i.e.,  0=0), 
the  associated  ionospheric  range  errors  are: 


1  -  fp/( 2 f 2 )  =  20.15  m 


(C.3-5) 


fpf  I cos6 | /f 2  =  2.82  x  I0'2m 


(C. 3-6) 


(C. 3-7) 


fp/(8f4)  =  1.01  x  l0'3m 
f2f2(sin20  +  2cos20)/(4f4)  =  9.87  x  10'6rn  (C.3-8) 

r  o 

Because  the  first-order  contribution  to  the  ionospheric  range 

error  is  four  orders  of  magnitude  greater  than  the  next  correc- 

2 

tion  term,  the  first-order  (1/f  )  contribution  to  the  ionospheric 
range  error  is  a  sufficient  approximation  for  a  radio  inter¬ 
ferometric  system  operating  in  the  gigahertz  frequencies. 

C . 4  CORRECTIONS  FOR  IONOSPHERIC  REFRACTIVITY  DELAYS 

To  correct  for  the  delay  of  electromagnetic  signals 
propagating  through  the  ionosphere,  upper  atmospheric  modeling 
or  dual- frequency  measurements  can  be  used.  Modeling  the  iono¬ 
sphere,  as  noted  in  the  previous  section,  requires  well-studied 
and  predictable  electron  behavior.  In  the  northern  mid-latitude 
regions ,  the  ionosphere  has  been  exhaustively  researched  and 
is  better  understood  than  the  polar  or  equatorial  areas.  Thus, 
in  these  well-studied  areas,  signal  delay  can  be  reduced  by  50 
to  70  percent  compared  with  uncorrected  measurements,  on  an 
rms  basis.  Alternatively,  since  the  delay  is  proportional  to 
the  inverse  square  of  the  signal  frequency,  two  frequencies, 
transmitted  coherently  and  sufficiently  far  apart,  can  also  be 
used  to  reduce  the  ionospheric  delay.  The  remainder  of  this 
section  addresses  the  uncertainty  associated  with  this  method. 

The  signal  path  difference,  p,  which  is  used  to  deter¬ 
mine  the  relative  receiver  distance,  is  shown  in  Fig.  C.4-1. 
Measurement  of  this  difference  yields 

m.  =  p  +  d  +  (K/f2)  +  n. 


(C.4-1  ) 


FIXED  SITE 
RECEIVER 

Figure  C.4-1  Interferometric  Signal  Path  Difference 

nr  =  signal  path  difference  measurement  using  channel 

d  =  tropospheric  and  mechanical  uncertainties, 
independent  of  frequency  and  receiver  noise 

K  =  parameter  proportional  to  the  difference  in  the 
total  electron  content  along  the  signal  paths 

f^  =  signal  frequency  of  channel  i 

n.  =  receiver  noise  in  channel  i. 


Tropospheric  and  mechanical  uncertainties  are  ignored  in  this 
discussion.  Measurement  in  two  channels  yields  an  equation 
for  the  true  signal  path  difference 

p  =  (m^  -  ni)  +  f^Km^-n^)  -  ]/( f2  -  f|)  (C.4-2) 

The  second  term  in  Eq.  C.4-2  stems  from  the  algebraic  elimina¬ 
tion  of  K,  whereas  the  first  term  is  the  single-channel  differ¬ 
ence  measurement  uncorrected  for  ionospheric  delays.  The 
receiver  noise,  n^ ,  is  an  unknown  by  definition,  so  the  esti¬ 
mate  of  the  true  path  difference  is 

p  =  +  f2(m^-m2 )/( f2- f^ )  (C.4-3) 

and  the  error  in  this  estimate  is  given  by 

6p  =  p  -  p  =  n^  +  )/(  fj"  )  (C.4-4) 

If  the  receiver  noise,  n^ ,  in  both  receivers'  channels 
have  equal  standard  deviations,  o^,  and  correlation  factor  p, 
then  the  standard  deviation  in  the  signal  path  difference 
error  is 

4^4  0  f 2.2, 1/2 

a  -  an(f2  +  fr2pflf2] 

=  an  {1  +  [^f1f2/(f|-ff  )  l2  ( 1-p  ) }  1/2  (C.4-5) 

The  first  term  in  the  braces  of  Eq.  C.4-5  represents  the  single 
channel  contribution  to  the  measurement  error.  The  second 
term  arises  from  the  ionospheric  calibration  process.  Fig¬ 
ure  C.4-2  is  a  plot  of  the  signal  path  difference  error  esti¬ 
mate,  Sp,  as  well  as  the  standard  deviation  in  the  signal  path 
difference  error  estimate  given  by  Eq .  C.4-5  for  the  case  in- 
\  >lvmg  uncorrelated  channels.  Using  the  current  GPS  trans- 
•’  i  1  t  a!  frequencies,  the  total  receiver  measurement  error  is 
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approximately  three  times  the  single-channel  receiver  noise. 
Although  the  ionospheric  delay  is  not  exactly  proportional  to 
the  wavelength  squared,  higher-order  terms  constitute  a  neg¬ 
ligible  source  of  error  (see  Section  C.l).  Ionospheric  signal 
delay  error  will  not  contribute  to  radio  interferometric  measure¬ 
ment  errors  until  the  inter-receiver  distances  are  large  enough 
to  cause  two  distinct  signal  paths  from  satellite  to  receivers 
(receiver  separations  are  greater  than  1000  km).  For  these 
very  long  distances,  dual - frequency  receivers  will  offer  cor¬ 
rections  to  the  signal  delay;  the  error  associated  with  the 
corrections  is  given  by  Eq .  C.4-5. 


APPENDIX  D 


THIRD-ORDER  MARKOV  PROCESSES  FOR 
ATMOSPHERIC  REFRACT IV ITY  MODELING 


D.l  INTRODUCTION 

Since  the  lower  atmosphere  is  not  static,  any  tropo¬ 
spheric  model  is  geographically  localized.  A  first  approxima¬ 
tion  to  tropospheric  modeling  requires  exponential  correlation 
of  site  atmospheric  conditions.  A  third-order  Markov  process 
was  selected  to  model  tropospheric  refractivity  effects  at 
varying  distances  from  a  weather  monitor  station.  This  choice 
was  made  to  ensure  distance  correlation  of  atmospheric  condi¬ 
tions  after  double-differentiation  of  inter-receiver  length; 
i.e.,  to  obtain  acceleration.  For  completeness,  a  full  deriva¬ 
tion  of  the  state-space  representation  for  a  third-order  Markov 
process  and  the  associated  covariance  matrices  is  presented. 
The  dynamics  matrix  is  derived  in  Section  D.2;  Sections  D.3 
and  D.4  define  the  spectral  density  and  initial  covariance 
matrices,  respectively. 


D.2  DERIVATION  OF  THE  DYNAMICS  MATRIX 

The  autocorrelation  fuction  for  a  critically  damped 
third-order  Markov  process  is  given  by  (Ref.  22) 

d>(T)  =  a2  e'PT[l  +  gt  +  Or)  2/3  ]  (D.  2-1  > 

where 


o  =  rms  error 

g  =  correlation  time 

r  =  time-shift  parameter. 


A  linear,  state-space  model  which  generates  the  random  process 
with  the  correlation  process  defined  by  Eq .  D.2-1  has  the 
following  form: 


where 


x  =  Fx  +  w 


(D.  2-2 ) 


x  =  state  vector 
F  =  dynamics  matrix 

w  =  vector  of  zero-mean,  white-noise  elements. 

To  derive  the  state-space  form  for  the  third-order 
Markov  process,  begin  by  defining  a  variable  x ^  as 

Xj_  =  e'pt  [1  +  pt  +  (pt)2/3]  (D.2-3) 

Taking  the  first  derivative  of  x^  (with  respect  to  t)  yields: 

xL  =  -pe‘px[pT/3  +  P2t2/3]  (D. 2-4) 

For  convenience,  a  variable  X2  can  be  defined  as 

x2  =  xL/p  (D.2-5) 

and  thus  , 

x2  =  -e"Px(PT/3  +  P2x2/3]  (D. 2-6) 

Equation  D.2-6  can  be  rewitten  in  terms  of  x,  by  replacing  the 

2  1 
t  term  with  an  equivalent  expression  obtained  from  Eq .  D.2-3. 

Thus  , 


=  -x 


+  e 


■PT 


(1  +  2  p  x / 3 ] 


( D  .  2  -  7  ) 


The  derivative  (again  with  respect  to  x)  of  x 3  yields: 
x 2  =  -x1  -  pe'PT[l/3  +  2pi/3]  (D.2-8) 

A  new  variable,  x^  ,  can  be  defined  as 

x3  =  *'c2//^  (D.2-9) 

so  that 

x3  =  -x1/p  -  e’^T [ 1/3  +  2pt/3]  (D.2-10) 

Rewriting  Eq.  D.2-10  explicitly  in  terms  of  x3  and  results 
in 

x^  =  -x^  -  2x3  +  2e  ^x/3  (D.2-11) 

where  use  has  been  made  of  Eq .  D.2-5  and  an  equivalent  expres¬ 
sion  for  the  x  term  was  obtained  by  rearranging  Eq .  D.2-7. 

The  time  derivative  of  is  given  as 

x^  =  -Xj  -  2x3  -  2pe  ^x/3  (D.2-12) 

or  equivalently, 

x3  =  ‘^xl  "  3^x2  ■  3gx3  (D.2-13) 

after  substituting  the  appropriate  expressions  for  x^ ,  x3,  and 
2e~®T/3  from  Eqs .  D.2-5,  D.2-9,  and  D.2-11,  respectively. 

Finally  arranging  Eqs.  D.2-5,  D.2-9,  and  D.2-13  in 
the  form  of  D.2-2  results  in 


D-3 


“ 

“  ■ 

-  - 

*1 

0  1  0 

X1 

x2 

=  p 

0  0  1 

x2 

1 - 

X  • 

l _ _ 

-1  -3  -3 

•  «■ 

.  X3. 

where  the  dynamics  matrix  is  given  by 


(D. 2-1-4 


F  =  P 


0 

0 

-1 


1  0 

0  1 

-3  -3 


(D.2-15) 


The  block  diagram  representation  for  this  third-order  Markov 
process  model  is  presented  in  Fig.  D.2-1. 
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F igure  D.2-1 


Third-Order  Markov  Process  Model 


D .  3 


THE  SPECTRAL  DENSITY  MATRIX 


The  white-noise  vector,  w,  from  Eq .  D.2-2  has  a  spec 
tral  density  matrix  Q  defined  by 

E[w(t)w(i)]  =  q  6(t-x)  (D.3-1) 

where  6(*)  is  the  Dirac  delta  function;  i.e.,  Q  is  diagonal. 
An  appropriate  expression  for  Q  can  be  written  as  (Ref.  22) 

0  0  0 

0  0  0 

0  0  16 


D . 4  INITIALIZATION  OF  THE  COVARIANCE  MATRIX 

The  covariance  of  the  state  vector,  denoted  by  P,  i: 
governed  by  the  differential  equation 

P  =  FP  +  PFT  +  Q  (D.4-1) 


where  superscript  T  indicates  the  matrix  transpose  operator. 
In  steady- state ,  P  =  0.  Hence,  defining  P  as 


P  =  cr‘ 


a  b  c 


d  e  f 


(D.4-2) 


in  steady  state,  Eq .  D.4-1  yields 


0 


d  e  f 

g  h  i 

-(a+3d+3g)  -(b+3e+3h)  -(c+3f+3i)_ 


+ 


b 

e 

h 


b  -(a+3b+3c) 

f  -  (d+3e+3f ) 

i  - (g+3h+3i ) J 


+  Q 


(D  • z 


-3) 


Solving  the  linear  equations  of  D.4-3  simultaneously  results 
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